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ABSTRACT 

Most, if not all, disk galaxies have a thin (classical) disk and a thick disk. In most models 
thick disks are thought to be a necessary consequence of the disk formation and/or evolution 
of the galaxy. We present the results of a study of the thick disk properties in a sample of 
carefully selected edge-on galaxies with types ranging from T = 3 to T = 8. We fitted one- 
dimensional luminosity profiles with physically motivated functions - the solutions of two stellar 
and one gaseous isothermal coupled disks in equilibrium - which are likely to yield more accurate 
results than other functions used in previous studies. The images used for the fits come from the 
Spitzer Survey of Stellar Structure in Galaxies (S 4 G). We found that thick disks are on average 
more massive than previously reported, mostly due to the selected fitting function. Typically, 
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the thin and the thick disk have similar masses. We also found that thick disks do not flare 
significantly within the observed range in galactocentric radii and that the ratio of thick to thin 
disk scaleheights is higher for galaxies of earlier types. 

Our results tend to favor an in situ origin for most of the stars in the thick disk. In addition 
the thick disk may contain a significant amount of stars coming from satellites accreted after the 
initial build-up of the galaxy and an extra fraction of stars coming from the secular heating of 
the thin disk by its own overdensities. 

Assigning thick disk light to the thin disk component may lead to an underestimate of the 
overall stellar mass in galaxies, because of different mass to light ratios in the two disk compo- 
nents. On the basis of our new results, we estimate that disk stellar masses are between 10% 
and 50% higher than previously thought and we suggest that thick disks are a reservoir of "local 
missing baryons" . 

Subject headings: galaxies: photometry — galaxies: spiral — galaxies: structure 

1. Introduction 

Thick disks in lenticular and spiral galaxies arc defined as disk-like components with a scaleheight larger 
than that of the thin or 'canonical' disk. They are typically detected as exponential excesses of light a few 
thin disk scaleheights above the galaxy mid-plane in edge-on galaxies. They have been first detected by 
Tsikoudi (1979) and recognized as a distinct galaxy structural component by Burstein (1979). A thick disk 
component was soon found in the Milky Way (Gilmore & Reid 1983) and it was afterwards proved to be 
made of old and metal-poor stars when compared to the thin disk stellar population (Reid & Majewski 1993; 
Chiba & Beers 2000). Later studies found thick disks to be nearly ubiquitous (Yoachim & Dalcanton 2006; 
Comeron et al. 2011a) and to be systematically older than their thin counterparts (Yoachim & Dalcanton 
2008a). 

The origin of thick disks is still a matter of debate, with four formation mechanisms being discussed (see 
the introduction of Yoachim & Dalcanton 2008a and references therein) . The first two possibilities consider 
two different mechanisms to dynamically heat an originally thin disk, meaning that its vertical stellar velocity 
dispersion is increased. In the first possibility this dynamical heating could be a consequence of internal 
evolution due to gravitational interaction with thin disk overdensities like giant molecular clouds or spiral 
arms (Villumscn 1985; Hanninen & Flynn 2002; Schdnrich & Binney 2009; Bournaud et al. 2009). In the 
second the heating is considered to be due to galaxy-galaxy or dark matter subhalo-galaxy interactions and 
mergers (Quinn et al. 1993; Robin et al. 1996; Walker et al. 1996; Velazquez & White 1999; Chen et al. 2001; 
Benson et al. 2004; Hayashi & Chiba 2006; Kazantzidis ct al. 2008; Read et al. 2008; Villalobos & Helmi 
2008; Qu et al. 2011). The third possibility is that the thick disk is a consequence of in situ star formation 
(Brook et al. 2004; Robertson et al. 2006; Elmegreen & Elmegreen 2006; Brook et al. 2007; Richard et 
al. 2010) or star formation with a high initial velocity dispersion in very massive star clusters (Kroupa 2002). 
Possibilities one and three are related if thick disk star formation occurs in gas that has a high turbulent 
velocity dispersion from its own clump stirring. The fourth possibility is that the thick disk is formed by 
accretion of stars from disrupted small satellite galaxies after the initial build-up of the galaxy (Statler 1988; 
Gilmore et al. 2002; Abadi ct al. 2003; Navarro et al. 2004; Martin et al. 2004; Read et al. 2008). 

Another topic of debate is that of "missing" or "lost" baryons (Persic & Salucci 1992). The fraction 
of baryons detected at high redshift (see, e.g., WMAP predictions by Spergel et al. 2003), does apparently 
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Fig. 1. — Reproduction of Figure 12 in van der Kruit & Searle (1981) using S 4 G data. The plot displays 
luminosity profiles of NGC 4565, going from a galactocentric distance R = outward using a radial bin 
width of 46". The data of the four quadrants have been averaged in order to produce the plot. 

not match the present-day fraction of baryons. Part of the missing baryons are expected to be found in 
hot gas in the warm-hot intergalactic medium (WHIM; see, e.g., Fukugita et al. 1998, simulations from Cen 
& Ostriker and WHIM emission detection/constraining by Hickox & Markevitch 2007 and Zappacosta et 
al. 2010). Additionally, dark matter haloes in which galaxies lie are also found to be very baryon deficient 
compared to the cosmic baryon fraction (McGaugh 2008). This apparent lack of baryons inside the galaxy 
potential well is what has been called the "local missing baryon" problem. Part of these baryons have 
probably escaped due to winds caused by supernova feedback, but the local baryon budget is far from being 
accounted for (McGaugh 2008 and references therein). 

In this Paper we make use of Spitzer Survey of Stellar Structure in Galaxies (S 4 G; Sheth et al. 2010) mid- 
IR imaging of edge-on disk galaxies to assess their stellar vertical mass distribution. The main observational 
goal of S 4 G is obtaining IRAC (Fazio et al. 2004) 3.6/xm and 4.5/Ltm-band images for 2331 galaxies with 
a radial velocity V TB dio < 3000 km s -1 . These passbands are not strongly affected by dust (Draine & Lee 
1984; Sheth et al. 2010) and are not as sensitive to star-forming regions as optical wavelengths, making 
them especially well suited for tracing the underlying stellar population of galaxies. In addition, having 
a typical surface brightness sensitivity of /iAB ~ 27 mag arcsec -2 in the 3.6/im band the S 4 G images have 
an unprecedented depth for such a large number of galaxies. This high sensitivity allows us to reproduce 
luminosity profiles down to /iAB ~ 26 mag arcsec -2 at 3.6/im (Figure[l} with a much better spatial resolution 
than presented in previous papers (e.g., Burstein 1979; van der Kruit & Searle 1981). In order to study these 
edge-on galaxies, we made fits of vertical luminosity profiles in the 3.6/im band. The fitted function results 
from coupling two luminous (stellar) isothermal disks in equilibrium. We also considered a third component 
made of cold gas in some cases. We tested the effect that the inclusion of a dark matter halo would have in 
our fits and found it to cause biases smaller than the error introduced by choosing the mass-to-light ratios 
of the thin and thick disk. In addition, we tested the effect of fitting the luminosity profiles of non perfectly 
edge-on galaxies and found it to introduce a small bias which goes opposite to that introduced by a dark 
matter halo. 

This Paper is structured as follows. We present the selected sample in Section 2. In Section 3 we 
describe the luminosity profile fitting we adopted and describe its potential caveats. The results of our fits 
are presented in Section 4 and are discussed in Section 5. The conclusions and a brief summary of the Paper 
are presented in Section 6. 
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Fig. 2. — 3.6/im-band S 4 G image of NGC 0522, one of the galaxies in our sample. NGC 0522 is similar to 
the Milky Way both in morphological type and in size. The vertical lines indicate the limits of the bins for 
which luminosity profiles have been produced, from left to right, at galactocentric distances of —0.8 r 25 < R < 
-0.5 r 25 , -0.5 r 25 < R < -0.2 r 25 , -0.2 r 25 < R < 0.2 r 25 , 0.2 r 25 < R < 0.5 r 25 and 0.5 r 25 < R < 0.8 r 25 . In 
order to avoid the influence of the bulge we have ignored the central (—0.2 r 2 5 < R < 0.2 r 2 s) bin throughout 
the Paper, in all galaxies. 



2. Sample selection 

The sample studied in this Paper was selected from the 817-galaxy S 4 G subsample for which reduced 
data were available at the time (December 2010; Salo et al. 2011). Ellipticity profiles were constructed 
during the S 4 G data processing from the 3.6/xm images by running iraf's ellipse with the galaxy center 
position fixed. We selected all those disk galaxies with types — 3 < T < 8 whose maximum disk ellipticity is 
e > 0.8. Galaxies of morphological type 8 < T < 9 (the so-called Magellanic galaxies) were rejected because 
their structure is generally ill-defined. The resulting 61 selected galaxies were visually inspected in order 
to detect any sign of spiral structure or (pseudo)rings, which would indicate a disk which is not edge-on. 
The images used for this rejection process were those from S 4 G as well as some from the Sloan Digital Sky 
Survey Data Release 7 (Abazajian et al. 2009) and the Hubble Space Telescope Legacy Archive (hosted at 
|http: / /hla.stsci.eduj ). As a result, fourteen galaxies were removed from the sample. We removed a final 
galaxy from the sample because of the presence of a very bright star close to the galactic disk affecting it at 
most galactocentric radii. 

The final sample contains 46 galaxies with types 3.0 < T < 7.5 at a median distance of D = 24.7 Mpc 
(using the NASA/IPAC Extragalactic Database mean value of redshift-independent distances except for 
IC 1553, NGC 5470 and PGC 012349 for which we estimated the distance using HyperLEDA's heliocentric 
radial velocity from radio measurement and using a Hubble constant Hq = 70kms~ 1 ). The sample lacks 
early-type galaxies for two reasons. First, a large fraction of the S 4 G galaxies are late-types (see Figure 1 of 
Sheth et al. 2010). Second, the ellipticity criterion used for selecting the sample biases against galaxies with 
large bulges or stellar haloes. 

As a comparison, one other recent study on thick disks in late-type galaxies with a statistically significant 
sample is that from Yoachim & Dalcanton (2006), who study 34 edge-on galaxies of types 5 < T < 9 with a 
median distance of D = 74.3 Mpc. 

An S 4 G 3.6/im-band image of one of the selected galaxies, NGC 0522, is shown in Figure O 
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3. Fitting procedure 

3.1. Choice of the fitting function 

Traditionally, single-component edge-on disks have been fitted by functions proportional to sech 2 (z/z ), 
where z is a scaleheight. This function was derived on theoretical grounds by Spitzcr (1942) and van der 
Kruit & Searle (1981) by solving the equations of equilibrium of a single isothermal sheet. In this context 
isothermal means that the stars in the disk are assumed to have a constant vertical dispersion velocity, 
((^z) 2 )' j at all z. This function can be approximated as an exponential with scaleheight z C xp = zq/2 at large 
vertical distance, which justifies the exponential fits used by several authors (Morrison et al. 1997; Comeron 
ct al. 2011a). Later on, due to discrepancies between the real luminosity profiles and the sech 2 (z/2:o) function 
at low z, van der Kruit (1988) proposed the rather ad hoc function sech 2 ^™(riz/2zo), as a generalization of 
his previously proposed fit. Recently Banerjee & Jog (2007) have suggested that these discrepancies are due 
to the gravitational interaction of the thin disk with the gas disk. 

The functions described previously are reasonable approximations for the behavior of a single disk, 
and superpositions of them have been used to describe edge-on galaxies with both a thin and a thick 
disk component (Shaw & Gilmore 1989; Shaw & Gilmore 1990; de Grijs & van der Kruit 1996; Morrison et 
al. 1997; Naslund & Jorsater 1997; Wu et al. 2002; Pohlen et al. 2004a; Yoachim & Dalcanton 2006; Comeron 
ct al. 2011a). However, Morrison et al. (1997) make a warning about the limitations of one-dimensional fits 
arising from degeneracies. While these superpositions are handy because they are analytic, they are not as 
physically justified as in the case of a single-component disk, as they ignore the gravitational interaction 
between the different disk components. This is why, for this Paper, we choose to integrate the equations of 
equilibrium for a set of gravitationally coupled isothermal stellar and gas disks. 



3.2. Creation of synthetic luminosity profiles 

We assumed that disks are relaxed structures whose stars or gas clumps behave like particles of an 
isothermal fluid in equilibrium. The assumption of equilibrium was recently proven to be reasonable for the 
Milky Way (Sanchez-Salcedo et al. 2011). We assume each disk (stellar thin disk, stellar thick disk and gas 
disk) has its own internal velocity dispersion. Under these conditions the equations which the disks follow 
are: 

d 2 pi pi ( . . dif D M\ 1 /dpi\ 2 

-47rG(p t + Pt + Ps) + — +— -j- , (1) 



where t refers to the properties of the thin disk, T to those of the thick disk, g to those of the gas disk, 
the subindex i can be either t, T or <?, p stands for the mass density, and ((f z ) 2 ) 5 for the vertical velocity 
dispersion of the particles - stars or gas clumps - as described by Narayan & Jog (2002). di ^° M is a 
term which describes dark matter effects. From now on we will define the vertical velocity dispersion as 
°~\ = (( w z) 2 ) 5 - To solve this set of three second-order differential equations we used the Newmark-/? method 
with j3 — 0.25 and 7 = 0.5, where (3 and 7 are internal parameters of the algorithm which have been set to 
make it unconditionally stable (Newmark 1959). 

We set the mid-plane to be at z = and then we integrated the equations from z = to a large z = Zf , 
where Zf causes that p(z = 0) » p(z = Zf), p being the sum of the densities of all the disks. For practical 
purposes Zf can be considered as infinity. In order to solve these three coupled second order equations six 
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boundary conditions are required. The selected ones were p\(z — 0) (pio thereafter) and ^-| z= o- The first 

da 

dz 



set of derivative boundary conditions is naturally ^P-\ z= o = as the mid-plane is by definition the location 



of the maximum density in the disk. 

We solved our differential equations for a grid of models with different central density ratios (ratios of 
initial conditions; pto/pw) and different vertical velocity dispersion ratios, or/ov We made integrations for 
150 values of PTo/Ptcb equally spaced from 0.015 to 2.25, and for 150 values of or/o^, equally spaced from 
f .f to 16.0. The limiting values of the grids were selected a posteriori, after testing grids of different sizes 
and verifying that the selected one was covering the whole parameter space of the studied sample. 

Two cases were studied, namely one with / Z Z =0 °° p g dz = and one other with f z ~£° p g dz = 0.2 / Z Z =0 °° Pt dz. 
The first case represents a disk with no gas and the second one represents a disk with an average abundance 
of gas (the gas mass fraction at Solar radius is 0.17 according to Banerjee & Jog 2007). We refer to these 
cases as 'without-gas' and 'with-gas', respectively. The gas velocity dispersion was fixed to be o~ g = (l/3)<Jt, 
in rough agreement with the values found for the disks of the Milky Way, of a g = 8kms _1 for the local Hi 
gas (Spitzer 1978), a g — Skins" 1 for the local Hn gas (Stark 1984; Clemens 1985), and a thin disk stellar 
velocity dispersion a t = 18kms _1 (Lewis & Freeman 1989; Narayan & Jog 2002; Banerjee & Jog 2007). 

We set the dark matter halo term in Equation 1 to be dJ ^" M = 0. This will be further justified in 
Section 3.6.2. 

We transformed the resulting synthetic stellar mass density profiles into synthetic luminosity profiles 
by assigning a different mass-to-light ratio (T) for the thin and the thick disk stellar components. In order 
to calculate T we used the Bruzual & Chariot (2003) spectral synthesis models with the set of stellar 
evolutionary tracks from the 'Padova 1994' library (Alongi et al. 1993; Bressan et al. 1993; Fagotto et 
al. 1994a; Fagotto et al. 1994b) with a Salpeter Initial Mass Function (IMF; Salpeter 1955). Once selected a 
given Tx/Tt, we made the synthetic luminosity profiles and scaled them to have a central luminosity equal 
to one and to p ^~—q^ =0.1, where the units of z are arbitrary. We considered three different star formation 
histories (SFHs) calculated for the solar neighborhood. These three SFHs are shown in Figure [3] and are 
described in the next subsection. 

In summary, the equations need six boundary conditions and a scaling factor for the absolute value 
of the intensity. The intensity scaling was done by setting the mid-plane intensity equal to 1 in both the 
models and the observations. Then the six boundary conditions consisted of three differentials at z = 0, 
i.e., -A| s= q = for each component, two normalizations, which are the thick-to-thin disk density ratio and 
velocity dispersion ratio, and one other normalization for the vertical pixel scale, which was taken to be 
£t*=2°°l=0.1. 

p(z=0) 



3.3. Adopted star formation histories 

3.3.1. Pilyugin & Edmunds star formation history 

Pilyugin & Edmunds (1996) used the following function to fit the stellar population of the Solar neigh- 
borhood, which we assumed to be representative of a typical thin disk population: 

a (i/ r to P ) exp(-t/r t0 p) for t < T top 



1 exp(-f/T SF R) for t > T top 
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Fig. 3. — SFHs inferred for the Solar neighborhood, from Pilyugin & Edmunds (1996; red), Nykytyuk & 
Mishenina (2006; green) and Just & Jahreifi (2010; purple). The horizontal axis denotes time with the origin 
13 Gyr ago. The vertical axis shows the star formation rate with an arbitrary scaling. Solid curves denote 
the thin disk SFH, and dashed curves denote the thick disk one. 



where a = exp(l — Ttop/TSFri) is a normalization constant and where Tt op and tsfr are free parameters which 
define the timescale of a first burst of star formation and the timescale of the subsequent star formation 
decay. We used the parameters obtained from their best fit, namely T top = 1 Gyr and tsfr = 8 Gyr for a 
13 Gyr old thin disk. Approximating the results from their chemical evolution model (Figure 2 of Pilyugin & 
Edmunds 1996) we set thin disk stars formed at time Gyr < t < 2.75 Gyr to have a metallicity of Z — 0.004, 
those formed at 2.75 Gyr < t < 7 Gyr to have Z = 0.008 and those formed at 7 Gyr < t < 13 Gyr to have 
Z = 0.02. As the authors of this study make no statement about the thick disk population, we assumed 
it to be formed quickly at the beginning of the galaxy history - in practice as a result of a single burst 
of star formation 13 Gyr ago (t = 0) and with a low metallicity (Z — 0.0004). As a result we found that 
T-r/Tt = 2.4. The calculated T accounts only for the stellar mass (stars and stellar remnants) and ignores 
the gas expelled from stars in the form of stellar winds, supernovae and planetary nebulae. This is not a 
problem since part of the gas was assumed to be expelled from the galaxy by the supernovae winds and the 
remaining gas was taken into account in the modeled gas component. We used the same IMF as Pilyugin & 
Edmunds (1996), namely that described by Salpeter (1955). 

As the star formation history is probably different for every galaxy, we have calculated how Tx/T t would 
vary with tsfr- We have found that Tx/T t = 2.0 for tsfr = 5 Gyr and Tx/T t = 2.5 for tsfr = 10 Gyr. 



Nykytyuk & Mishenina (2006) used the same SFH - i/j(t) -, IMF, and Galaxy age as Pilyugin & Edmunds 
(1996) for the thin disk. For the thick disk they used as similar function with tsfr = 5 Gyr. Approximating 
the results from their chemical evolution model (Figures 5 and 8 of Nykytyuk & Mishenina 2006) we set 
thin disk stars formed at time OGyr < t < 2.75 Gyr to have a metallicity of Z = 0.004, those formed at 
2.75 Gyr <t<7 Gyr to have Z = 0.008 and those formed at 7 Gyr < t < 13 Gyr to have Z = 0.02. For the 
thick disk we set the stars formed at time Gyr < t < 0.75 Gyr to have a metallicity of Z = 0.0004, those 
formed at 0.75 Gyr < t < 1.25 Gyr to have Z = 0.004 and those formed at 1.25 Gyr < t < 13 Gyr to have 
Z = 0.008. For this SFH we calculated T T /T t = 1.2. 



3.3.2. Nykytyuk & Mishenina star formation history 
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3.3.3. Just & Jahreift star formation history 

In the model of Just & Jahreifi (2010), the Milky Way is assumed to have an age of 12 Gyr and a SFH 
as follows: 

(t + t )tl 

where to = 5.6 Gyr, ii = 8.2 Gyr and t n = 9.9 Gyr. Approximating the results from their chemical 
evolution model (Figure 16 of Just & Jahreifi 2010) we set thin disk stars formed at time Gyr < t < 0.75 Gyr 
to have a metallicity of Z = 0.004, those formed at 0.75 Gyr < t < 7.5 Gyr to have Z = 0.008 and those 
formed at 7.5 Gyr < t < 12 Gyr to have Z = 0.02. Just & Jahreifi (2010) adopted the thick disk to have 
an age of 12 Gyr and here we have set it to have a metallicity Z = 0.0004. For this SFH we calculated 
Tx/Tt = 2.4. This value is only an approximation as Just & Jahreifi used a Scalo (1986) IMF when fitting 
the stellar population of the solar neighborhood and using a Salpeter (1955) IMF as we have done introduces 
some internal incoherence in the calculations. 



3.3.4. Which SFH should we use? 

We can summarize the three previous subsections by stating that there are basically two kinds of 
SFH (see Figure 2). The first class of models implies that the stars in the thick disk were created nearly 
instantaneously early in the galaxy history like in our implementation of a thick disk in the Pilyugin & 
Edmunds (1996) model and in Just & Jahreifi (2010). The two models yield a very similar result, namely 
that Tx/Tt = 2.4. The second class of models implies that the thick disk is made of stars created over a 
longer time, like what is proposed by Nykytyuk & Mishenina (2006) yielding Tx/T t = 1.2. Choosing one 
of these two kinds of SFHs would imply some assumptions on the formation mechanism of the thick disk. 
Moreover, different galaxies are likely to have different star formation histories which will yield a different 
T T /T t . 

A rough estimate of Tx/T t can be made from the data collected in Yoachim & Dalcanton (2006). The 
authors used their B — R measurements for the thin and the thick disks in order to estimate (Tx/T t ) fl using 
some recipes included in Bell & de Jong (2001), a paper which also includes recipes to convert B — R to 
(Tx/Tt) K . As the X-band is relatively close to the 3.6/zm-band, (Tx/Tt)^ can be used as an approximation 
for the (Tx/Tt) in this paper. For the galaxies in Yoachim & Dalcanton (2006), (Tx/T t ) K = 1.6 with a 
dispersion of 0.3. Using their data we found that there is a certain trend of (Tx/T t ) K with the maximum 
circular velocity of the galaxy, v c , with (Tx/T t ) K = 2.1 ± 0.2 — (0.005 ± 0.002)w c . The linear correlation 
coefficient is r — 0.5. This relationship means that more massive galaxies tend to have disks with a lower 
(Tx/T t ) K . 

In order to avoid a Tx/T t choice, and also due to the large scatter in the Tx/T t derived from Yoachim 
& Dalcanton (2006), we decided to study two cases, namely Tx/Tt = 1.2 and Tx/Tt = 2.4, which are the 
results of different fitted Milky Way SFHs. The Tx/Tt of most galaxies is probably found between these 
values, as discussed in the previous paragraph. As we also studied two gas disk concentrations, we produced 
four grids of 150 by 150 luminosity profile models. The Tx/T t ratio affects the thick and thin disk stellar 
mass ratios by nearly a constant factor. In what follows, for many figures, we only show the Tx/Tt = 1.2 
results and then mention the appropriate factors for scaling these results to the T T /T t = 2.4 fits. 
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3.3.5. Limitations and caveats of the selected fitting function 

We assumed the thin and the thick disk to be made of isothermal sheets with two distinct and unique 
vertical velocity dispersions. Although this is true in first approximation, a small but significant increase 
of the vertical velocity dispersion has been found for Milky Way thick disk stars at large vertical distance 
above the mid-plane (Moni Bidin et al. 2010). 

The line of sight integration has not been taken into account when producing the synthetic luminosity 
profiles. This is a valid approximation as long as the radial scalelengths of all the disks are similar and the 
scale heights have no radial dependence. The first assumption allows the ratios of mid-plane densities of 
the different disks to be constant. In the sample studied by Yoachim & Dalcanton (2006) the thick to thin 
disk scalelength ratio is on average (/it/Th) ~ 1.2, which justifies our assumption. The second assumption 
prevents the integration of different scale heights on a line of sight that covers a range of radii. Our results 
here are consistent with this assumption, as we see no significant change in the thick disk scaleheight with 
radius. 

The third caveat is that we neglected the effect of the bulge. This can be justified because our sample 
is made of late-type galaxies and because the fits were made at R > 0.2 r25, where the bulge is relatively 
unimportant. 

We also ignored the effects of dark matter and not perfectly edge-on galaxy geometries, which will be 
justified in Sections 3.6.2 and 3.6.3. 

3.4. Compilation of observed luminosity profiles 

Using S 4 G 3.6/im-band images, we produced four luminosity profiles for each galaxy: at each side 
of the bulge along the long disk axis for projected galactocentric distances 0.2 r^h < \R\ < 0.5 T25 and 
0.5r25 < \R\ < 0.8r25, as shown in Figure [2] (r25 from HyperLEDA; Paturel et al. 2003). The profiles were 
prepared by adding the counts above and below the mid-plane in order to get one unique profile for each 
bin and averaging for each z over the non-masked pixels. We used manually refined masks coming from the 
S 4 G Pipeline 2 (Sheth et al. 2010). 

3.5. Comparison between observed and synthetic luminosity profiles 

The observed luminosity profiles were scaled in luminosity and in height above the mid-plane, z, using 
the same procedure as for the synthetic profiles. We ignored the mid-plane pixel in our fits to avoid the 
region of highest extinction. Note that at 3.6/im, and especially at R > 0.2 r25 where we study vertical 
profiles, there is very little dust extinction. However, to allow for a small amount of extinction, we scaled 
the luminosity profile outside of the mid-plane by various factors, from 0.80 to 1.00 in steps of 0.01. These 
20 different scalings for the luminosity profiles introduce a third dimension to the parameters of the fit. The 
synthetic profiles were then convolved with a Gaussian kernel with a Full Width at Half Maximum (FWHM) 
of 2.2", which is the FWHM of a synthetic 3.6/im-band S 4 G Point Spread Function (PSF), made by stacking 
several hundreds non-saturated stars (S 4 G "super-PSF"; Sheth et al. 2010). Each pixel in S 4 G has a size of 
0.75". 

The code which compared the observed and synthetic luminosity profiles was made of two nested loops 
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Fig. 4. — Flow diagram showing the procedure used to compare the real and the synthetic luminosity profiles, 
as described in Section 3.5. 



which arc illustrated in the flow diagram shown in Figure |4j 

The external loop compared the brightest section of the observed luminosity profiles, with different 
assumed dust extinctions, with a Gaussian convolved synthetic profile. The brightest section was defined 
to have a given dynamical range, Am, with an initial Am = 3 mag arcsec -2 . We term the faintest fitted 
magnitude as The best fitted profile - that with the lowest sqrt\ 2 in magnitudes - was considered to be 
that with the 'correct' central dust extinction. If the best fit had y^x 2 " < 0.1 mag and if (i\ < 26 mag arcsec -2 , 
we made an extra iteration of the loop with Am = Am + 0.5 mag arcsec -2 . If one of these conditions was 
not fulfilled, we considered as the best fit that obtained for Am = Am — 0.5 mag arcsec -2 . The limit 
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Fig. 5. — Fits to the surface brightness profiles of NGC 0522 for the four fitted bins for the case Tr/Tt = 1.2 
and with-gas. The solid curve represents the observed luminosity profile, and the dashed curve the best fit. 
The dotted curves indicate the contributions of the thin and the thick disk. The dash-dotted vertical lines 
indicate the limits of the range in vertical distance above the mid-plane used for the fit. 



jiil = 26 mag arcsec 2 was selected in order to avoid any noise effect at low surface brightness regions. 

The inner loop was in charge of detecting the thickness of an eventual mid-plane dust layer. The 
synthetic and the real luminosity profiles were compared for z > 1 pixel in order to avoid comparing regions 
which are presumably affected by dust. We selected as a best fit the synthetic profile for which the \px^ was 
minimum. Then we compared again the real profiles with the synthetic grid, but for z > 1.2 pixels (we made 
linear interpolation between pixels in order to use fractional pixels). If the best fitting synthetic profile was 
the same as when comparing for z > 1 pixel, we took this fit as a good result; if not, we considered that 
at z = 1 pixel the dust effect is still important and thus repeated the procedure comparing the best fits for 
z > 1.2 pixels and for z > 1.4 pixels. We continued this procedure with steps of Az = 0.2 pixels until the 
best fit for z — z a was the same as for z — z a + Az. 

The statistical uncertainties in the fits are negligible compared to those introduced by considering a 
given Tx/Tt. 

An example of the fitted profiles is presented in Figure [5] 
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Fig. 6. — Plot of several of the synthetic surface brightness profiles in the without-gas case. The panel on 
left shows profiles computed with Tx/T t = 1.2 and those at the right with Tx/T t = 2.4. The horizontal 
axis shows the vertical distance above the mid-plane in the arbitrary units we have used for our scaling. The 
vertical axis shows the decrement of surface brightness in magnitudes relative to the mid-plane. Red, green 
and violet curves indicate profiles with pTo/pto equal to 0.09, 0.77 and 2.25, respectively. Short-dashed, 
long-dashed and solid curves denote profiles with or/fit equal to 1.76, 2.47 and 3.62, respectively. The 
horizontal black lines indicate, from top to bottom, Am equal to 4.0, 4.5 and 5.0 magnitudes. 




Fig. 7. — Plot comparing the synthetic surfaced brightness profiles in the without-gas case with those in the 
with-gas case. The profiles have Tx/T t = 1.2 and pto/pw = 0.77, like the green profiles in the left panel in 
Figure HI The green curves denote the without-gas case and the violet ones the with-gas case. 



3.6. Reliability of the fits 

3.6.1. Degeneracies 

We tested our fits in order to detect possible degeneracies and to test the conditions in which our fitting 
method can be applied. 

Figure [5] shows the synthetic profiles obtained for a selection of pTo/pto and i^t/oV For pTo/pto = 0.77, 
a central density typical of what can be found in our fits, the profiles with a different ratio of velocity 
dispersions are easily distinguishable because the inflection point which indicates the transition from the 
thin to thick disk-dominated area can be detected even when the fit is produced over a low Am (where 
Am is the range of magnitudes over which the fit is produced). The profiles are much more degenerate 
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Fig. 8. — Square root of the mean square difference between a single-disk profile and our grid of synthetic 
two-disk profiles. The vertical axis displays the ratio of the mean vertical dispersion velocity between the 
thin and the thick disk, ax/crt, and the horizontal axis displays the ratio between the thin and thick disk 
mid- plane density, pxo/pto- The selected Am is given in the bottom- left corner of each panel. All plots are 
for Tx/T t = 1.2 and with-gas. For Tx/T t — 2.4, the distribution is very similar, but shifted by a factor 
two to the right. The with-gas and without-gas distributions are nearly identical. 



for pto/pw — 2.25, the case with a higher central thick disk density in our grid, as the luminosity is thick 
disk-dominated at every z, hiding the inflection point. At the other side of the studied range of central 
density ratios, the profiles made for thin disk-dominated galaxies - Pto/pw = 0.09 - show a clear inflection 
point, but only if Am is large enough. In the case of Tx/Tt = 2.4, a galaxy with pxo/pto = 0.09 shows a 
profile very similar to that of a galaxy with pxo/pto — 2.25 unless Am becomes large enough to unveil the 
inflection point. 

Figure [7] shows how adding a gas disk in the model of a galaxy affects the synthetic luminosity profiles 
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for a model with an intermediate pTo/pto- The effect is small, and causes the model thin disk scaleheight to 
become smaller. 

To quantify the possibility of a two-disk galaxy to be identified as a single-disk one, we estimated 
the mean square difference between our sets of synthetic luminosity profiles and an isothermal single-disk 
profile in equilibrium for different total magnitude ranges in the fit, Am (Figure [8]). We show only the case 
Yt/1\ =1.2 with-gas; for Tx/T t = 2.4, the distribution is very similar and shifted a factor two to the 
right. The with-gas and without-gas distributions are nearly the same. We have considered a profile to 
be compatible with a single-disk luminosity distribution when yx 2 " < 0.2 mag arcsec -2 . As indicated by 
Figure [6j we have found that for high values of pto / Pto the profiles are degenerate and a two-disk profile 
cannot be distinguished from a single-disk profile. When Am is low, this degeneracy also affects profiles 
with a low pto/ Pto- In order to avoid degeneracies, we have decided to include in our study only galaxies for 
which the luminosity profiles have been successful over a dynamical range of Am > 4.5 mag arcsec -2 . This 
restriction, combined with that of fitting until the faintest level of p\ = 26 mag arcsec -2 , ensures that only 
profiles with a central brightness, fi(z — 0) < 21.5 mag arcsec -2 , are considered. Two galaxies, even though 
they have p(z = 0) < 21.5 mag arcsec -2 , are excluded from the study, because the fit could only be done 
over Am < 4.5 mag arcsec -2 due to their profile not being compatible with that of two coupled isothermal 
stellar disks in equilibrium. 



3.6.2. Effect of ignoring a dark matter halo 

We ignored the effect that a dark matter halo would have in our fits. Including them would be prob- 
lematic due to our ignorance of the properties of the dark matter haloes, with possibilities ranging from 
maximum disks (for which no significant amount of dark matter amount is needed to explain the shape 
of the rotation curve within the optical radius) to clearly submaximal disks, going through studies which 
indicate that disks with large circular velocities are maximum and those with lower circular velocities are 
submaximal (see discussion in Bosnia 2004). Examples of studies which have found at least disks in high- 
luminosity galaxies to be maximum are those by Salucci & Persic (1999) and Palunas & Williams (2000). 
Moni Bidin et al. (2010) have found no significant dynamical effect of dark matter in our Galaxy at Solar 
radius, implying that its disk is likely to be maximal. On the other hand, studies over statistically significant 
samples which find indications of an opposite result include those by Courteau & Rix (1999) and Pizagno et 
al. (2005). 

We produced a "worst-case" test on how ignoring the presence of a dark matter halo affected our fits. 
As the main results of this paper relates to the ratio of column mass densities between the thin and the thick 
disk, Ex/St, our tests were directed to know how this parameter varies with the inclusion of a dark matter 
halo. To do so we have made several tests on NGC 5470. We choose this galaxy because it is the faintest 
one in our sample for which it has been possible to do fits with Am > 4.5, and it thus is the studied galaxy 
which is more likely to have a submaximal disk. We have used the Narayan & Jog (2002) formalism, with a 
pseudo-isothermal dark matter halo (van Albada et al. 1985) term as follows: 

dKDM <aX r i-K~ ,fI C * +^( 3 4-l) -tan , (4) 



dz r 2 \ r 2 r 2 (R 2 + r 2 ) r \ r 2 J \R C 

where u ma x is the extrapolated rotation speed at infinity, R c is the core radius of the dark matter halo and 
r = y/ p 2 + z 2 is the 3D radius or distance to the galaxy center. 

As the inclusion of a dark matter halo implies the need for a line-of-sight integration, we needed to know 



- 15 - 



Table 1: Column density ratios (Ex/St) obtained from fitting NGC 5470 luminosity profiles including the 
effect of a dark matter halo. 
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Note. - The fit made with w max = 0(a) corresponds to that described in Section 3.5. The fit made with 
timax = 0(b) corresponds to a fit made with no vertical pixel scaling, such as in the case of fits with dark 
matter, implying an extra fitting parameter (at and ctt instead of at/ar)- 

the scalelength of the disk. NGC 5470 has a disk truncation (downbending) at R <~ 40", which corresponds 
to R <~ 3kpc. The bins with 0.5 r25 < \R\ < 0.8 r25 - those which are more likely to be affected by the dark 
matter halo - are located in the outer, steeper disk. Using van der Kruit & Searle (1981) approximations 
for measuring the scalelength of edge-on disks when R is much larger than the scalelength, we found that 
the scalelength of the outer disk is h <~ 700 pc. We also found the galaxy to have no significant emission in 
3.6^m for R > 6.2 kpc, a distance we have set as the end of the stellar disks. In addition, we assumed that 
the mid-plane stellar mass density of NGC 5470 at p = 0.65r25 is p t o + Pto = 0.1 Mq pc~ 3 , which is a value 
on the order of the mass density in the Solar neighborhood. In order to reduce the computation time we 
tested without-gas cases only. The inclusion of a dark matter halo forces us to use the thin and the thick 
disk vertical velocity dispersions as free parameters (at and ctt) instead of their ratio, as we have been doing 
until now. 

As we aimed to use realistic dark matter halo parameters, we searched for galaxies similar to NGC 5470 
in de Blok et al. (2008), where the authors deduce halo properties from Hi rotation curves. We found that 
NGC 7793 has a luminosity and a rotation speed similar to that of NGC 5470, so we used the values of 
R c corresponding to this galaxy found in their Table 3 (R c — 1.93 kpc and R c = 3.52 kpc) for testing the 
effect of dark matter. We considered cases with a maximum circular velocity v c = v max <~ 120 km s -1 and 
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\v c = f m ax ~ 60 kms~ . In addition we made tests with the halo of a much larger galaxy, NGC 0925, which 
has R c = 16.86 pc. We also considered several disk scalelengths around the value h ~ 700 pc, as we have 
found our results to be dependent on this parameter. The results are presented in Table [T] 

We found that, in general for a submaximal disk, by neglecting the halo we introduce a relatively small 
bias, namely we overestimate the ratio of Ex/S t by 10 — 30% However, for some peculiar haloes, the effect 
may be larger, up to ~ 40 — 50%. Considering that NGC 5470 is the faintest galaxy in the sample and that 
brighter galaxies are probably less affected by dark matter (closer to maximum) we conclude that ignoring 
the effect of a dark matter halo introduces an bias smaller than that caused by choosing a given Tx/T t 
(which has an effect on the order of 50% for the range of "reasonable" Tx/T t we have been studying), and 
that it is thus justified to ignore its influence. 

3.6.3. Effect of the line-of -sight integration in not 'perfectly edge-on galaxies 

Even though we have been quite restrictive at selecting galaxies to be included in our sample, some 
galaxies may be far from an ideal edge-on geometry. Thus, we tested the effect of line-of-sight integration in 
galaxies that are not perfectly edge-on. 

The test has been done for a "typical" galaxy in our sample, namely a galaxy with a radius r 2 5 = 14.3 kpc, 
a truncation radius somewhat larger than r25 , rtrunc — 20 kpc, a thin and thick disk scalelength h = 2 kpc and 
a thick disk scaleheight zt = 700 pc. In addition we set the distance to the test galaxy to be D — 25 Mpc. All 
the galaxy properties are scalable, and the distance to the galaxy has only been included to account for the 
effects derived from the PSF. We assumed T T /T t = 1.2 and selected the galaxy to have PTo/Pto = 0.50 and 
o~T/o~t = 2.19 which are also typical values obtained in the fits to the galaxies in our sample (see Section 4). 

The result of the test is shown in Figure [9] in which luminosity profiles have been fitted for different 
inclination angles, i, at the galactocentric radii, R, corresponding to the middle of the 0.2 r 25 < \R\ < 0.5 r 25 
and the 0.5 r 25 < \R\ < 0.8 r 25 bins and for Am = 5.0magarcsec -2 . The effect of lowering the inclination 
is to increase the fitted pto / Pto and to reduce the fitted ctt /ft , which translates into the fitted thin and 
the thick disks being less differentiated until to a certain i for which the sum of a thin and a thick disk is 
indistinguishable from a single-disk profile (the square root of the mean square difference between a single- 
disk profile and the fitted synthetic luminosity profile is y^X 2 " < 0.2 mag arcsec -2 ). The change in the fitted 
Pto/ Pto and or/ft happens in such a way that Ex/S t remains roughly constant. The small changes in 
Ex/St go in the direction of underestimating it (an effect on the order of ~ 10%), thus partly counteracting 
the effect of a dark matter halo in the case of a submaximal disk. The fitted thick disk scaleheight (zt) 
remains roughly constant when decreasing i, but the fitted thin disk scaleheight (z t ) increases significantly. 

We have repeated the test for a variety of reasonable galaxy parameters. In general two disks can be 
distinguished down to 80° < i < 85°. The fitted Ex/E t is in general in agreement with that of an edge-on 
geometry down to at least i = 86°, although it tends to be underestimated by a factor up to ~ 20%. For 
lower i values Ex/E t can be heavily underestimated, especially when measured at high R for galaxies with 
low r trunc (e.g. rtrunc = ^25) and for all R in those galaxies with low pTo/pto and/or low or/crt- 
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Fig. 9. — Plots representing the effect of fitting the luminosity profiles of a galaxy (with the parameters 
appearing in the bottom-left corner) which is not perfectly edge-on. From bottom to top, the solid lines 
represent how the luminosity profile would appear for decreasing inclination angles starting at i = 90° and 
with a step of At = 1°. The dotted lines represent the fitted thin and thick disk components. The table 
in the top-right corner shows parameters derived from the fits for different inclination angles. R stands for 
the galactocentric distance at which we are making the fit and (x 2 ) 1 ^ 2 is the square root of the mean square 
difference between a single-disk profile and the fitted synthetic luminosity profile in units of mag arcsec -2 . 
The values of R we selected correspond to the middle of the 0.2 r25 < \R\ < 0.5 r25 bin (left panel) and 
the middle of the 0.5 r25 < \R\ < 0.8 bin (right panel). The luminosity profiles have been plotted for 
decreasing inclination until y^ 2 = 0.2 mag arcsec -2 , which is value for which we considered a two stellar 
disk structure to be indistinguishable from a single-disk structure. 



3.6.4. Testing the selected Tx/T t 

A way to test whether the two Tx/T t chosen are two good limiting cases is to calculate the Tx/T t at a 
wavelength other than 3.6^m and compare the results with those obtained using S 4 G 3.6^m images. If the 
adopted SFH, and thus the Tx/T t , was accurate, then Sx/^t should be similar in both cases. 

We produced the test using SDSS r-band imaging for NGC 5470. We choose this galaxy because it 
has a thin and symmetric dust lane, thus facilitating the analysis. Using the SFHS in Section 3.3 and 
the Bruzual & Chariot (2003) spectral synthesis models we calculated that when Tx/T t = 1.2 in 3.6/xm, 
(Tx/T t ) r = 1.1 and that when Tx/T t = 2.4, (TT/T t ) r = 1.3. As our fitting code was designed to 
deal with low extinction levels and the r-band image of NGC 5470 is strongly affected by a mid-plane 
dust lane, we manually set the code to ignore the pixels for which z < 4". The fact that the dust has a 
large effect for more than half of the range in z where the thin disk luminosity dominates illustrates the 
importance of infrared imaging for this kind of study. We have been able to produce satisfactory fits for 
galactocentric bins except for that at 0.5 r25 < R < 0.8 T25- We found that on average when the 3.6/im 
T T /T t = 1.2, (S T /St) r = 1.42(E T /St) and that when T T /T t = 2.4, (E T /St) r - 0.78(E T /St), with a 
scatter compatible with the fitting uncertainties. This proves that for NGC 5470 the actual Tx/T t (that 
which makes (Ex/Et^ = (Ex/E t )) is found somewhere between the two limiting values used in this Paper. 
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3.6.5. The case of NGC 4565 

We found an argument which, although maybe circumstantial, seems to prove that using the solution 
of equations of equilibrium at fitting the luminosity profiles gives a more accurate result than a sum of 
somewhat ad hoc analytical functions. One of the most studied edge-on galaxies, NGC 4565, has been 
decomposed several times. Jensen & Thuan (1982) and Naslund & Jorsater (1997) needed three disks to fit 
it (or two disks plus a "corona" or halo component), but we only needed two to fit the luminosity profiles 
of this galaxy down to 26 mag arcsec -2 , which is similar to the depth of previous studies. The simpler fit 
certainly advocates for the use of non-analytic physically based functions as done in this study. Alternatively, 
the need for a third disk may be introduced by some peculiar dust geometry which we avoided by studying 
the galaxy in the infrared. 

4. Results 

Of the 46 galaxies selected to be fitted, 14 have a mid-plane brightness too low to be fitted down to 
Am = 4.5 mag arcsec -2 at any of the vertical bins that we have studied. In addition, two more galaxies, 
although bright enough, could not be fitted down to Am = 4.5 mag arcsec -2 because their profiles differed 
significantly from any of our synthetic models. Thus, we obtained good fits for only 30 galaxies in our sample. 
The parameters of the fits with Am > 4.5 mag arcsec -2 are presented in Tables [5] and [3] 

Most of the fits are compatible with a double stellar disk. However, some of the fits for ESO 548-063, 
IC 5052, NGC 1827, PGC 013646 and UGC 10297 are compatible with a single-disk structure (see discussion 
in Section 3.6.1). All these fits have Am = 4.5 mag arcsec -2 , so they have been calculated over a relatively 
small dynamic range, which may explain the non-detection of a thick disk. Alternatively, these galaxies may 
not be inclined enough (i ;$ 85) for allowing the detection of a two stellar disk structure. The results of fits 
compatible with a single disk have not been included in the plots of this Section. 

Two galaxies, ESO 079-003 and NGC 4013, have some bins with lower fitted magnitudes of fi\ < 
24 mag arcsec -2 . An inspection to their luminosity profiles and the images of these galaxies shows the 
presence of a third extended component, maybe a third disk or an especially bright stellar halo. NGC 3628 
also appears to have an extra extended component but faint enough to allow fitting the luminosity profiles 
down to levels fainter than /ij < 24 mag arcsec -2 . NGC 4013 is a peculiar galaxy in the sense it is known 
to have a very prominent Hi warp (Bottema et al. 1987) and a 'giant tidal stream' (Martinez-Delgado et 
al. 2009) and it has been studied in detail in a follow-up letter (Comeron et al. 2011b). For the other galaxies, 
our set of synthetic fitting functions appears to describe reasonably well all the disk components. 

Although for the same luminosity profile fits with different Tx/T t and gas content usually have the same 
Am, in a few cases fits including a gas disk seem to fit slightly better the luminosity profiles. For Tt /T t = 1.2, 
in the without-gas case the luminosity profiles were fitted over (Am) = 5.12 mag arcsec -2 on average, while 
when including a gas disk they were fitted over (Am) = 5. 19mag arcsec -2 . In the case of Tx/T t = 2.4 the 
average dynamical range of the fits is (Am) = 5.06 mag arcsec -2 and (Am) = 5.16 mag arcsec -2 , respectively. 
These calculations have been performed over the luminosity profiles for which the fit was successful for the 
two values Yx/T t and the two gas disk fractions used. As discussed by Bancrjcc & Jog (2007), this is 
probably due to the fact that the inclusion of an invisible gas disk helps to fit better the inner parts of the 
stellar luminosity profile under the assumption of equilibrium. 

We have found that for most cases the dust extinction in the mid-plane is small or negligible, as can be 



Fig. 10. — Left panel: central thick to thin disk mass density ratio as a function of galaxy morphological 
type. Right panel: thick to thin stellar column mass density ratio as a function of galaxy morphological type. 
Circles denote values fitted at 0.2 r 25 < \R\ < 0.5 r 2 5 and triangles values fitted at 0.5 r 25 < \R\ < 0.8 r 2 §. 
Both plots are for fits made using Tx/Tt = 1.2, with-gas. For Tx/Tt = 2.4, points are shifted on average 
by a factor 1.9 upwards in both plots. Plots for fits in the without-gas case are similar to those with-gas. 

seen in columns 19 to 22 of Table [2j 

There is no variation in the range of possible disk density ratios as a function of galaxy type, T, as can 
be seen in the left panel of Figure 1101 Moreover, we find that the area populated by points corresponding 
to 0.2 < \R\ < 0.5 (circles) is very similar to those corresponding to 0.5 r 2 $ < \R\ < 0.8 r 2 $ (trian- 
gles) showing that, as previously assumed when not taking into account the line of sight integration, the 
scalelengths of the thin and the thick disk are similar. 

The shape of the Ex/E t distribution is very similar to that of the one of pro/PtOj as can be seen in 
the right panel of Figure [TU] We find that log (Ex/E t ) scales linearly with log (pto/Pw), with a correlation 
factor of 0.58. Thus, the ratio of velocity dispersions, ax/ct, only accounts for a relatively small scatter. 

From the data plotted in FigureHUlwe find that for Tx/T t = 1.2 only ~ 60% of the data points indicate 
a higher column mass density in the thin disk than in the thick disk (Ex/E t < 1; ~ 70% of the points if an 
20% effect on Ex/E 4 is introduced by the halo as described in Section 3.6.2). This fraction goes down to 
~ 30% for Tx/Tt = 2.4 (~ 35% of the points if a 20% halo bias is considered). In some cases, the column 
mass density of the thick disk appears to be two or even three times larger than that of the thin disk. As a 
comparison, the baryonic column mass of the thick disk in the Milky Way has been estimated to be around 
~ 20% (Gilmore & Reid 1983; Chen et al. 2001; Robin et al. 2003; Juric et al. 2008). If we consider the gas 
fraction of the disk to be one order of magnitude less than that of the stars (Banerjee & Jog 2007), the Milky 
Way would appear as a very thin-disk dominated galaxy if Tx/Tt = 1.2, and as an outlier if Tx/Tt = 2.4. 
On the other hand, one other Milky Way study, by Fuhrmann (2008), considers the mass of Milky Way's 
thick disk to be comparable with that of the thin disk, which would make it a 'typical' galaxy in our sample. 
If we compare our results with estimates of thick disk masses in external galaxies we see that Yoachim & 
Dalcanton (2006) found that for several galaxies (seven out of 34 or ~ 20%) the stellar mass of the thick 
disk is greater than that of the thin disk, but for these galaxies, they assume a large gas fraction in the thin 
disk, including the far outer parts beyond the optical disk, making its total baryonic mass at least two times 
greater than that of the thick disk for all the galaxies in their sample. 

The ratio of vertical velocity dispersions, ot /<Jt, shows some correlation with morphological type, with 
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Fig. 11. — Left panel: central thick to thin velocity dispersion ratio in the z direction as a function of galaxy 
type. Right panel: ratio of thin and thick disk scaleheight as a function of the galaxy type. Circles stand 
for values fitted at 0.2 r 2 5 < \R\ < 0.5 and triangles stand for values fitted at 0.5 r 2 5 < \R\ < 0.8 r 2 5- 
The cross symbols denote averages and their errors in bins T = 1 wide. Both plots are for fits made using 
YT/Yt = 1-2) with-gas. For Tx/T t = 2.4, points are shifted on average by a factor 1.05 up in the velocity 
dispersion plot. Plots for fits in the without-gas case are similar to those with-gas. 
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Fig. 12. — Scaleheight for the bin 0.5 r 25 < \R\ < 0.8 r 25 , divided by that for the bin 0.2 r 25 < \R\ < 0.5 r 25 . 
Circles stand for the thin disk, and triangles for the thick disk. Plots for fits using Tx/T t = 2.4 and 
without-gas are similar to that presented here. 



the disks in earlier- type galaxies more kinematically differentiated (left panel in Figure ITT]) . We find that 
the degree of differentiation is not linked to the galaxy brightness or to v c . 

The right panel of Figure [TTJ shows the ratio of the scaleheights between the thin and the thick disk, 
Zt/z\. The scaleheight of each disk was measured by finding the distance between those points where the 
density of the disk is pi(z) = e _4 pio and Pi(z) = e _5 pio, where the subindex i stands for the thin or the thick 
disk. In the selected regime, the slope of each component is close to exponential, thus our scalelengths can 
be compared to those of exponential fits and to half of the scaleheight which is found in sech 2 (z/2:o) fits. 
The distribution of disk scaleheights is naturally very similar to that of vertical velocity dispersions, as, at 
first approximation, the velocity dispersion scales with the scaleheight of the disk. As zt is very sensitive to 
the inclination angle, i, part of the scatter in the plot may be caused by the not exactly edge-on orientation 
of the galaxies. 
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Fig. 13. — Ratio of the thick to thin disk stellar mass as a function of the circular velocity. The plot in the 
left panel is produced for Tx/T t = 1.2, and that in the right panel with Tx/T t = 2.4. The circle indicates 
a very asymmetrical galaxy, and crosses indicate galaxies with an X-shape bulge (see text). Both plots are 
in the with-gas case and should be compared to Fig. 22 of Yoachim & Dalcanton (2006). 



Disks, both thin and thick, do not systematically flare in a significant way, as seen in Figure [12l The 
only points at which the ratio Zi(0.5r25 < \R\ < 0.8r 2 5)/;zi(0.2r25 < \R\ < 0.5 r 2 s) is significantly higher 
than unity correspond to IC 2135 and NGC 4013. In the case of IC 2135, the flare appears linked to a warp. 
What we see as a flare in NGC 4013 could be a consequence of a bad fit, as NGC 4013 is one of the two 
galaxies in our sample which is likely to need more than two disks to be well fitted (see discussion at the 
start of this Section and Comeron et al. 2011b). 



5. Discussion 

5.1. Comparison of thick to thin disk mass ratios with Yoachim &; Dalcanton 

A recent study of a statistically significant number of thin/thick disk decompositions was presented 
by Yoachim & Dalcanton (2006). Here, we will compare our relative thick disk masses with their results. 
As mentioned in Section 4, they found that ~ 20% of galaxies have a thick disk stellar mass larger than 
that of the thin disk. We have not directly measured the disk masses, so in order to obtain an estimate of 
this quantity which could be compared with their results we have estimated the disk masses by using the 
following formula: 

where the subindcx 6=1 — 4 refers to the different bins in galactocentric distance for which the fits have 
been done and fi c = /ii — Am is the mid-plane surface brightness for a given bin. The results are presented in 
Figure [T3l where we plot Mr/M t against the maximum gas circular velocity as obtained from HyperLEDA 
(Paturel et al. 2003). We see that for T T /T t = 1.2, 11 out of 27 galaxies (~ 40%) have a thick to thin disk 
stellar mass ratio larger than unity and that this number goes up to 18 (~ 70%) when Tx/Tt = 2.4 (the 
sample of galaxies for which fits were fitted is 30, but three were dropped from this analysis as the fits in all 
the bins are compatible with a single disk mass distribution). 

The two plots in Figure [T3] present a hint of a trend with rapidly rotating galaxies having lower relative 
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thick disk mass (lower M^/Mt) than slowly rotating ones. This trend appears in a much sharper way in 
the Fig. 22 of Yoachim & Dalcanton (2006). In order to infer whether the large spread in our Figure [T51 
may be influenced by sample selection biases, we looked for particularities in our sample galaxies. We found 
one galaxy - IC 1553 - which is significantly distorted and asymmetrical. We also found four galaxies with 
X-shaped bulges, ESO 079-003, ESO 443-042, NGC 3628 and NGC 4013. Once these galaxies are removed 
from Fig. [T3l the trend discovered by Yoachim & Dalcanton (2006) appears clearly. The outlier behavior 
of the galaxies with X-shaped bulges remains even after considering only the bins at galactocentric radius 
0.5 r25 < \R\ < 0.8 T25, implying that their position in the plot is not due to their bulge being extended and 
affecting some of our fits. In addition, the four galaxies with an X-shaped bulges are known to be warped 
(Sanchez-Saavedra et al. 1990; Sanchez-Saavedra et al. 2002; also seen in S 4 G images). However, the warp 
alone does not cause these galaxies to be outliers, as other galaxies in the sample, such as NGC 0522 and 
NGC 4565, are also warped. 

Why do thick disks in galaxies which host an X-shaped bulge not fit the 'normal' behavior of Mx/M t ? 
Yoachim & Dalcanton (2006) selected their sample to be made of bulgeless galaxies, so galaxies with an 
X-shaped bulge were not included. However, our sample has not been selected to be bulgeless and contains 
several galaxies with a significant bulge which fit in the general trend. For example NGC 0522 and NGC 4565 
have boxy bulges (but not X-shaped) and fall into the main relationship. Two of the X-shaped bulges 
reside in the two galaxies for which they are bins with luminosity profiles that we could not fit down to 
Hi > 24.0magarcsec _1 , so their position in the plot could be a consequence of a bad fit, and in the case 
of NGC 4013 it is a galaxy which is accreting external material. NGC 3628 has been fitted down to a low 
Hi, but has obviously been recently perturbed as indicated by its tidal tail (Kormendy & Bahcall 1974). 
ESO 443-042 has no noticeable features except for its bulge and warp. So, in at least two out of four cases, 
the X-shaped bulge may have been caused by bar creation and buckling triggered by a minor merger event, 
as suggested by Mihos et al. (1995). Alternatively, the X-shape could be the manifestation of a spiral-like 
perturbation going down to the center of the galaxy, triggered by a perpendicular interaction (Figs. 6 and 7 
in Elmegreen et al. 1995). Other authors, such as Binney & Petrou (1985) have also suggested mergers as a 
cause for boxy bulges. An alternative, and more widely followed, possibility is that the X-shaped bulges are 
edge-on views of strong bars (Athanassoula & Misiriotis 2002; Martnez-Valpuesta et al. 2006) and do not 
involve any merging. They would then fall in the same category as the boxy/peanut bulges, the difference 
being either due to the viewing angle of the bar (e.g. Bureau et al. 2006), or to the strength of the bar (see 
Athanassoula 2008, for a review). Thus weaker bars seen edge-on would have a boxy shape and stronger ones 
an X shape, while stronger bars could look boxy-like if viewed sufficiently close to end-on. However, finding 
a link between the mechanism responsible for the X-shaped bulge and the atypical behavior of the thick disk 
is not straightforward, although it may be argued that the gravitational influence of a recent minor merger 
or a strong bar may cause us to observe transient non-equilibrium states for which our equilibrium fitting 
functions may not hold. 

If we ignore galaxies hosting an X-shaped bulge, our fraction of galaxies in which the thick disk contains 
more mass than the thin disk is 8 out of 23 (~ 35%) for Tx/T t = 1.2, and 15 out of 23 (~ 65%) for 
T^T/ft = 2.4. This remains in sharp contrast with the 20% found by Yoachim & Dalcanton (2006). One 
possibility is that this difference is due to the choice of the fitting function. Yoachim & Dalcanton's (2006) 
results for each galaxy come from doing the median of the results of six fits with different mid-plane maskings, 
different smoothings to account for the seeing and different combinations of sech(z/zo) and sech 2 (z/zo) 
accounting for the thin and the thick disk. Of these six fitting models, three make use of a sum of two 
sech 2 (z / zq) functions and are thus likely to dominate the median of the six models. 
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Fig. 14. — Comparison of the fit in the bin — 0.5r25 < R < — 0.2r25 for NGC 0522, for the equations of two 
stellar and one gaseous isothermal coupled disks in equilibrium (left) and the addition of two sech 2 (z/zrj) 
functions (right). Solid curves represent the observed luminosity profile, and the dashed curves the best fit. 
The dotted lines indicate the contributions of the thin and the thick disk. The dash-dotted vertical lines 
indicate the limits of the range in vertical distance used for the fit. 
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Fig. 15. — Comparison of the column densities obtained from fitting the luminosity profiles by integrating 
the equations of two stellar and one gaseous isothermal coupled disks in equilibrium with those from fitting 
the addition of two sech 2 (z/zo) functions. The dotted line indicates the best linear fit which crosses the 
origin. Circles stand for values fitted for 0.2 r 25 < \R\ < 0.5 r 25 , and triangles stand for values fitted for 
0.5 r 25 < \R\ < 0.8 r 25 . 



We tested whether fitting the sum of two sech 2 (z/zo) functions gives more weight to the thin disk than 
the functions we used. We produced sech 2 (z/zrj) fits over our luminosity profiles using the same procedure 
than when fitting the equilibrium solutions, except for the fact that, for simplicity, we have not used a PSF 
convolution and we have not taken into account the effect of dust extinction which we show to be very small 
in most cases (see Table We found that the shape of the thick disk close to the mid-plane appears rounded 
in the case of a sech 2 (z/zo) fit, but it is significantly more peaked when making a fit using the equations of 
equilibrium (see Figure [Ml for an example). 

From a linear regression crossing the origin (Figure fT5|). we found that (Ex/St) s = 0.59±0.03 (Et/SiOu, 
where s denotes the fits using two sech 2 (z/zo) functions, and h stands for the fits using the equations of 
equilibrium. This implies that, if we assume our approach to be accurate, and Yoachim & Dalcanton (2006) 
mass measurements to be dominated by the models with a sum of two sech 2 (z/zo) functions, they might 
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have underestimated Mx/Mt by a factor of ~ 1.7. When multiplying by 1.7 the Mx/M t values in Fig. 22 
of Yoachim & Dalcanton (2006), we find that 12 out 34 galaxies (~ 35%) would have Mr/M t > 1, a value 
close to that we obtain with Tx/Tt = 1.2. 

The gradient in Tx/Tt with v c discussed in Section 3.3.4 is likely to have some effect in the trend seen 
in Figure [T3l As Tx/Tt is larger for galaxies with a small v c , it would cause the trend to be strengthened. 

For further comparison with the results of Yoachim & Dalcanton (2006) we should decompose galaxies 
from their sample with our methodology using S 4 G imaging. There is a four galaxy overlap between the 
S 4 G and their sample but, unfortunately, at the time of the sample selection of this paper (December 2010) 
those galaxies have not gone through the automated mosaicking/cleaning process yet. We plan to study this 
overlap in detail in a follow-up paper. 

In summary, we find that our thick to thin disk mass ratios (Mx/Mt) are higher than those found in 
literature due to the large fraction of mass attributed to the thick disk at low z when fitting luminosity 
profiles with the equations of coupled isothermal disks in equilibrium. These results have been obtained 
from data which can reproduce Yoachim & Dalcanton's (2006) results when fitting a sum of two sech 2 (z/zo) 
functions. As our model is based in physically motivated functions and we are using images less affected by 
dust than in previous studies, our results are likely to be the most accurate so far for such a large sample. 
We will now explore the implications of our results, in particular higher thick disk masses. 

5.2. Thick disk formation mechanisms 

The main result of this Paper is that the total stellar mass of thick disks is similar to that of thin disks 
and that in some cases, the thick disk mass can be up to 2 - 3 times larger than that in the thin disk. These 
surprisingly massive thick disks, together with the fact that we have found no significant flares, put some 
constraints on their formation mechanisms. 

The first formation mechanism to be discussed is that of a dynamical heating of an originally thin disk 
by disk overdensities. Secular heating due to giant molecular clouds (GMCs) has sometimes been suggested 
as a reason of disk heating. However, this mechanism has been recognized as insufficient to explain the 
highest velocity dispersion stars (Villumsen 1985; even a fourfold increase in the number of GMCs would not 
be enough to explain the presence of the hottest populations of stars according, to Hanninen & Flynn 2002). 
This mechanism is thus not likely to create very massive thick disks. Alternatively, the disk may have been 
mostly heated in a short interval of time at the beginning of the history of the galaxy, an epoch in which 
there is a large population of clumpy disks whose clumps formed through disk gravitational instabilities and 
had a mass of a few 10 8 Mq (Elmegreen & Elmegreen 2005). These massive clumps have a stellar scattering 
power much larger than present-day GMCs. Bournaud et al. (2009) found that this mechanism plus further 
secular evolution could cause a thick to thin disk mass ratio Mx/M t ~ 0.3 in simulations designed to produce 
galaxies with Milky Way mass. Thus, Bournaud's et al. (2009) Milky Way-like simulated galaxies, although 
a little bit short in thick disk mass, would be compatible with the v c — Mx/Mt relationship for Tx/Tt = 1.2. 
However, as these simulations have only been run for Milky Way-like galaxies, the effect caused in a less 
massive galaxy in unknown. Thick disks created through star scattering and radial mixing by spiral arms 
has been discussed and modeled by Schonrich & Binney (2009) who also succeed to reproduce a reasonable 
thick disk (although again light when put in the v c — Mx/Mt relationship due to the fact that they fitted 
the Milky Way which we have shown to have often been measured to be very thin disk-dominated). 
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Another possibility is that the thick disk has been created by the heating of a thin disk by satellite 
galaxies or dark matter subhaloes. This mechanism is more efficient for galaxies with orbits far from coplanar 
with the disk (Read et al. 2008). According to simulations, these interactions may heat the disk to form a 
thick disk, but they are thought to also cause a significant flare whose exact properties are dependent on the 
encounter parameters (Quinn et al. 1993; Walker et al. 1996; Kazantzidis et al. 2008; Villalobos & Helmi 
2008; Bournaud et al. 2009). Most of these authors suggest the flare to be at galactocentric radii small 
enough to be detectable in some cases. The exception is the work by Villalobos & Helmi (2008), where the 
flare would only affect the thick disk at very large galactocentric radii, which would need deeper photometry 
than ours to be detected. Kazantzidis et al. (2008) suggested that the flares in edge-on galaxies are what 
we interpret as truncations, as flares would be so substantial that the stars would be scattered over a large 
vertical range, causing the disk in those regions to have a very low surface brightness when seen edge-on. 
This would also imply that truncations are more often detected in edge-on galaxies, which is likely to be 
the case (see discussion in Pohlen et al. 2004b). We have found that strong flares are not common among 
our galaxies, implying that thick disks are unlikely to form through the heating of a disk by satellites. Even 
if the flares were so strong to appear as truncations, the heating through encounters would fail to explain 
the existence of thick disks in some galaxies, since thick disks are ubiquitous and truncations, although very 
common, do not appear in all galaxies (Kregel et al. 2002; van der Kruit 2007). Minor mergers indeed occur 
in disk galaxies, but they are likely to have a small effect at high redshift, where they are more frequent, 
due to the presence of a high gas fraction in the disk (Moster et al. 2010; previous simulations usually do 
not consider gas effects). The effect of gas for avoiding heating by interactions with satellite galaxies is still 
significant for galaxies with a gas fraction similar to that of present-day Milky Way, although it allows part 
of the thin disk to be heated into a flared thick disk (Qu et al. 2011). 

Some authors theorize that a thick disk may be a consequence of an in situ star formation. According 
to the simulations of Brook et al. (2004), the thick disk star formation is triggered by the accretion of 
gas-rich protogalactic fragments during hierarchical clustering at high redshift. The gas disk would still be 
dynamically hot at the moment of the thick disk creation, which would explain the difference in scaleheight 
between the thick disk which formed early and the thin disk formed from a cooled gas disk and from additional 
gas coming from cold flows. Brook et al. (2004) found that thick disk scalelengths are larger than those of 
thin disks, which is in contradiction with observations (Yoachim & Dalcanton 2006). Later models using 
disk-disk gas-rich mergers (Robertson et al. 2006; Brook et al. 2007; Richard et al. 2010) are able to produce 
reasonable thick disks from stars from the progenitor galaxies plus stars created during the merger, with 
the stars which are created later settling in a thin disk. However, these simulations do not include posterior 
accretion of cold gas, which will settle in the thin disk, increase its mass, and reduce the scaleheight of the 
thick disk. 

Another in situ formation model is that from Elmegreen & Elmcgrccn (2006), who study chain and spiral 
galaxies in the Hubble Space Telescope Ultra Deep Field (UDF) and find that the scaleheight of those galaxies 
is similar to that of present-day thick disks, with a mass column density of 4 to 40 Mq pc~ 2 . Recent studies 
(e.g., Moni Bidin et al. 2010) found the Milky Way column mass density in the solar neighborhood to be 
~ 60Mq pc~ 2 , thus, if we assume UDF thick disks to be progenitors of present-day thick disks, thick-to-thin 
disk column mass density fractions up to Ex/St ~ 2 would be compatible with the UDF thick disk mass 
estimates for Milky Way-like galaxies. However, interpreting these galaxies as direct progenitors of thick 
disks for which a thin disk has been added after accreting gas through cold flows is problematic, as discussed 
in Elmegreen & Elmegreen (2006) as the UDF thick disks are likely to shrink to lower scaleheights when 
further material is accreted, making them much thinner than observed present-day thick disks. Elmegreen 
& Elmegreen (2006) calculated that in order to create a Milky Way-like thick disk with zt = 875 pc, the 
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equivalent UDF thick disk should have zt ~ 3 kpc, which is far larger than the observed average UDF thick 
disk. However, as we found that thick disks are much more massive than previously thought, the fraction of 
added gas needed to generate a thin disk is lower than that estimated by Elmegreen & Elmegreen (2006), 
making the decrease in zt of the UDF thick disks lower than what they estimated. 

In a ACDM scenario, at least a fraction of the thick disk stars is expected to be accreted from satellite 
galaxies dragged into the disk plane by dynamical friction. This process would be effective for galaxies with 
an impact angle 9 < 20°, and larger angle encounters produce spherical- like distributions as well as a thick 
disk by thin disk heating, as discussed before (Read et al. 2008; but again these simulations do not include 
gas and thus the effect of disk heating may be overestimated). However, according to these authors, this 
mechanism applied to the Milky Way would produce a disk ~ 2 - 10 times less massive than observed, and 
thus needs to be combined with other thick disk formation mechanisms. Scannapieco et al. (2011) found that 
their simulated thick disks, which have been defined to be made of stars older than 9 Gyr, have "relatively 
high fractions" of non-accreted stars although the exact value is not specified. A strong indication that this 
mechanism is likely to be at the origin of a significant fraction of material in at least some galaxies is the 
discovery of a thick disk with a significant fraction of counterrotating material by Yoachim & Dalcanton 
(2008b). 

We now summarize our comments on the various proposed thick disk formations scenarios: 

• Heating of a thin disk by its own overdensities, in the form of, e.g., spiral arms, GMCs or giant clusters: 
This mechanism may have difficulties to produce massive enough thick disks. 

• Heating of a thin disk by satellite galaxies or dark matter subhaloes crossing it: Not likely to cause a 
significant effect as it would cause noticeable flares. Its effect is highly reduced when the disk contains 
a significant fraction of gas. 

• In situ thick disk formation: This formation mechanism had the problem of the shrinking of the thick 
disk scaleheight due to the accretion of external material after the thick disk formation, which would 
make thick disks indistinguishable from their thin counterparts. However, our results indicate that 
thick disks of many galaxies are much more massive than previously thought, lessening the impact of 
the newer material in the thick disk scaleheight. 

• Accretion of stars of satellite galaxies (tidal stripping): This mechanism occurs necessarily for some 
galaxies, as indicated by their kinematical or morphological signatures. However, simulations have 
shown that this effect can not explain all the mass found in thick disks. 

Thus, our results plus data from previous studies tend to favor an in situ origin for most of the stars in 
the thick disk, where the thin disk created afterwards is not massive enough to reduce zt by too much. In 
addition the thick disk may contain a significant amount of stars coming from satellites accreted after the 
initial build-up of the galaxy and an extra fraction of stars coming from the secular heating of the thin disk. 

5.3. Thick disks: lair of missing baryons? 

Thick disk-dominated galaxies may play some role in what has been called the "missing baryon problem" 
(Persic & Salucci 1992). WMAP cosmic microwave background observations have constrained the fraction 
of baryons in the Universe (Spergel et al. 2003), but most of this baryonic matter is apparently not found 



-27- 



in the "obvious" reservoirs which are stars and cold gas in galaxies (Bell et al. 2003). At a galactic scale, 
several studies indicate that observations are missing a large fraction of the baryons relative to the dark 
matter content in the Milky Way halo, and this problem appears to be even more severe for galaxies with 
smaller masses (Hoekstra et al. 2005; McGaugh 2007; McGaugh 2008; Bregman 2009). Rough counts by 
McGaugh (2008) for the Milky Way show that we are missing at least 10 11 Mq, which is comparable or 
larger than the mass of all the baryonic mass thought to be contained in the disk. 

Two mechanisms have been proposed to explain the "missing baryons" in galaxy dark matter haloes. 
The first is that there is a quantity hot gas in the halo which has either failed to cool down, or was blown 
away by supernova feedback. This hot gas in the halo has been measured to have a mass of ~ 4 x 10 s Mq 
(Bregman & Lloyd-Davies 2007) , which falls way short of solving the "local missing baryon" problem. The 
second possibility is that supernova feedback has blown away halo gas so efficiently that this material has 
been lost into the intergalactic medium. The amount of baryons lost in this way is difficult to estimate but 
it is also apparently not enough to explain the low fraction of baryons in galaxies (McGaugh 2010). Even 
though they remain insufficient, both explanations have the merit of explaining satisfactorily why lower-mass 
galaxies have lower baryon fractions, as supernova feedback is more effective in shallower potential wells. 

For external galaxies, the stellar masses inside the disks are calculated by assuming a given mass-to-light 
ratio. However, ignoring the presence of an unexpectedly high fraction of stars in the thick disk may have 
led to an underestimate of the value of T. Some numbers may be easily obtained from the data plotted in 
Figure IT3l Let us consider the more thin disk-dominated galaxies (Mr/Mt = 0.3) for the case Tx/T t = 1.2. 
If we were ignoring a thick disk, and assign a thin disk T for all the light emitted by the galaxy, we would be 
underestimating the stellar mass of the galaxy by ~ 10%. On the other hand, if we consider the most thick 
disk-dominated galaxy (Mr/Mt = 3.4), and Tx/T t = 2.4, we would underestimate the stellar disk mass by 
a factor of two. These numbers are calculated for the 3.6//m-band, which implies that the underestimation of 
the disk stellar mass would be different when measured at shorter wavelengths. We stress that these numbers 
only account for the baryons in stellar disks, and that if extended gas disks beyond the optical radius were 
included in the thin disk mass as done by Yoachim & Dalcanton (2006) Mr/Mt would go down. 

Thus, our results suggest that a least a fraction of the "local missing baryons" are located in the thick 
disk component of galaxies and that they have been "lost" until now due to the underestimation of Mr/Mt. 
This explanation is complementary, and not an alternative, to those which have already been suggested to 
solve the problem. 

6. Summary and conclusions 

The nature of the mechanisms creating the thick disks has been a topic of discussion for a long time. 
Basically four pictures have been suggested, each leaving some distinct observational signature. The first 
one is that the thick disk is the consequence of the thin disk have been dynamically heated by its own 
overdensities. The second possibility is that thin disks being dynamically heated by the disk crossing of 
satellite galaxies or dark matter subhaloes. The third possibility is in situ star formation during or shortly 
after the build-up of the galaxy. The fourth is that the thick disk could be made of stars which have been 
accreted from infalling satellites after the formation of the thin disk. 

In order to study the problem of thick disk formation we have used 3.6/zm images of edge-on galaxies 
from the S 4 G. The selection criterion used, maximum isophote ellipticity e > 0.8, has led to a sample of 
46 late- type galaxies with types ranging from T = 3 to T = 8. We have fitted the luminosity profiles of 
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these galaxies in radial bins with 0.2 r 2 5 < \R\ < 0.5 r 2 5 and 0.5 r 25 < \R\ < 0.8 r 2 5, using solutions of the 
equations of equilibrium for two stellar and one gaseous coupled isothermal disks. These functions have the 
advantage of being physically motivated. The largest uncertainty of the fit comes from the assumed ratio 
of mass-to-light ratios of the thick and the thin disk, Tx/T t . We studied cases with Tx/T t = 1.2 and 
Tx/Tt = 2.4 in order to account for the range of fitted star formation histories for the Solar neighborhood. 
The low value of Tx/T t arises when the stars in thick disks form over a long period of time and the high 
value arises when the thick disk stars formed entirely in a burst in the early Universe. The studied interval 
in T-r/Tt is in broad agreement with that derived from the study by Yoachim & Dalcanton (2006). 

We found that the stellar mass of the thick disk component ranges from one third to two times that of 
the thin disk if Tx/Tt = 1.2 and from two thirds to more than three times if Tx/Tt = 2.4. We thus found 
that thick disks are much more massive than generally assumed so far. In order to check the reasons of this 
difference, we compared our results in detail with those from Yoachim & Dalcanton (2006). We found that 
the choice of the fitting function is responsible for the high fraction of stellar light we detect in the thick disk. 
We also found that disks do not flare significantly within the range of galactocentric radii we have studied, 
and that the ratio of thick to thin disk scaleheights is higher for earlier-type galaxies in our sample. The 
lack of flares discards the kinematical heating of a thin disk by satellite galaxies or dark matter subhaloes 
as the main thick disk formation mechanism. 

As a result of these findings we suggest the following tentative formation mechanism for thick disks. 
First, a thick stellar disk is created during and soon after the build-up of the galaxy by the stirring action 
from massive disk clumps and the instabilities that form those clumps. A thin disk subsequently forms 
from gas which has not initially been spent in stars, and from gas that arrives later at a lower rate so that, 
when combined with the hot disk and bulge that are already present, it is relatively stable and cannot heat 
itself much. The newly accreted gas makes the thick disk scaleheight, Zt, shrink, and if this shrinkage is a 
relatively small amount, then the final thick disk remains a distinct component of the galaxy, with a younger 
thin disk inside of it. In the case of thick-disk dominated galaxies, the thick disk forms partly at high redshift 
by internal disk stirring, as above, and then increases its mass over a Hubble time by the heating of the thin 
disk by internal disk substructure, such as large star clusters and giant molecular clouds, and the accretion 
of satellite galaxies. 

As our results suggest that thick disks have masses comparable to those of thin disks and in addition 
they have a large mass-to-light ratio due to the absence of young stars, we suggest that ignoring them has 
led to an underestimate of the baryonic mass located in disk galaxies. We thus suggest that thick disks 
contain part of what has been called the "local missing baryons" . We have calculated this underestimate to 
be between 10% and 50% of the stellar mass of a disk, depending on the ratio of the thick versus thin disk 
masses and mass-to- light ratios, Mx/M t and Tx/T t . 
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Table 2. Galaxy and fit data (part one) 
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1.95 


2.07 






1.00 


0.96 






IC 5052 


7.1 


8.1 


141.4 


425 




79.8 


V 
V 
V 


V 
V 
V 


V 
V 
V 


V 
V 
V 


0.95 
1.65 
1.71 


1.08 
1.80 
1.88 




1.84 
2.00 
1.95 


2.00 

2.14 
2.10 






1.00 
1.00 
1.00 


0.96 
0.96 
0.96 




















V 


V 


V 


V 


0.29 


0.29 0.33 


0.33 


2.17 


2.26 


2.37 


2.10 


1.00 


1.00 


1.00 


1.00 


NGC 0522 


4.1 


35.0 


32.9 


144 


-20.69 


169.1 


V 
V 


V 
V 


V 
V 


V 
V 


0.25 
0.69 


0.27 0.32 
0.60 0.60 


0.32 
0.61 


2.07 
2.35 


2.14 
2.45 


2.21 
2.47 


2.02 
2.21 


1.00 
1.00 


1.00 
1.00 


1.00 
1.00 


1.00 
1.00 
















V 


V 


V 


V 


0.50 


0.48 0.59 


0.60 


2.17 


2.24 


2.35 


2.12 


1.00 


1.00 


1.00 


1.00 
















V 


V 








0.25 






2.14 








1.00 






NGC 0678 


3.0 


27.1 


77.4 


185 


-20.95 


169.0 


V 
V 
V 


V 
V 
V 


Pd 


Pd 

ps 




0.10 
0.48 
0.23 






2.39 
2.24 
2.43 








1.00 
1.00 
1.00 




















_L 


V 


V 


V 




0.46 0.50 






2.51 


2.47 






1.00 


1.00 1 




NGC 1351A 


4.5 


20.9 


132.1 


147 


-18.17 


89.6 


_L 
_L 
_L 


V 
V 
V 


V 
V 
V 


V 
V 
V 




0.46 0.51 
0.83 0.89 
0.86 0.91 






2.39 
2.65 
2.55 


2.37 
2.61 
2.51 






1.00 
1.00 
1.00 


1.0000 

1.00 ! 
1.00 


















_L 


V 


V 


_L 




0.56 0.80 






2.17 


1.87 






1.00 


1.00 




NGC 1495 


5.0 


17.3 


104.5 


131 


-18.62 


90.5 


_L 
_L 

_L 


V 
V 
V 


V 
V 
V 


_L 
_L 
_L 




0.57 0.81 
1.00 0.44 
1.02 0.46 






2.07 
2.28 
2.19 


1.79 
1.95 
1.87 






1.00 
1.00 
1.00 


1.00 
1.00 
1.00 


















V 


V 


V 


V 




0.16 1.23 






1.90 


2.07 






1.00 


1.00 




NGC 1827 


5.9 


11.6 


119.6 


194 


-18.68 


77.1 


V 
V 
V 


V 
V 
V 


V 
V 
V 


V 
V 
V 




0.12 1.29 
1.18 2.16 
0.21 2.25 






1.92 
1.97 
2.00 


2.00 
2.12 
2.07 






1.00 
0.97 
1.00 


1.00 
1.00 
1.00 




















V 


V 


V 




0.45 0.53 






1.82 


1.87 






1.00 


1.00 




NGC 3501 


5.9 


23.3 


28.0 


262 


-19.11 


133.6 


* 


V 
V 
V 


V 
V 
V 


V 
V 
V 




0.39 0.53 
0.84 0.96 
0.78 0.95 






1.73 
1.90 
1.82 


1.79 
1.95 
1.87 






1.00 
1.00 
1.00 


1.00 
1.00 
1.00 


















_L 


V 


V 


_L 




0.48 0.53 






2.68 


2.59 






1.00 


1.00 




NGC 3628 


3.1 


12.2 


102.6 


658 


-21.44 


211.4 


_L 
_L 

_L 


V 
V 
V 


V 
V 
V 


_L 
_L 
_L 




0.48 0.56 
0.86 0.93 
0.86 0.98 






2.57 
2.85 
2.74 


2.49 
2.72 
2.63 






0.98 
1.00 
0.98 


1.00 
1.00 
1.00 





Table 2 — Continued 



ID 


T 


13 


PA 


D 25 




v c 




Status 






Pto/ 












fo 










(Mpc) 


(°) 


(") 




(kms" 1 ) 












Pta 




















(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(91 


do) 


(11) 


(12) 


(13) 


(14) 


(15) 


(16) 


(17) 


fl8) 


(19) 


(20) 


(21) 


(22) 


(23) 
















_L 


V 


V 


V 




0.29 


0.30 


0.25 




2.85 


2.90 


3.00 




1.00 


1.00 


1.00 


In 4U.LO 


3 


18 6 


65 1 


294 


-19 38 


181 7 


_L 
_L 


V 

± 


V 
_L 


V 
_L 




0.29 


0.29 


0.30 




2.70 


2.76 


2.85 




1.00 


1.00 


1.00 
















± 


V 


V 


V 




0.54 


0.56 


0.50 




2.93 


2.95 


3.15 




1.00 


1.00 


1.00 
















V 


V 


V 


V 


0.48 


0.33 


0.32 


0.33 


2.43 


2.17 


2.19 


2.32 


1.00 


1.00 


1.00 


1.00 


wpp A^^n 


6 3 


19 5 


59 2 


137 


-19 87 


115 5 


V 
V 


V 
V 


V 
V 


V 
V 


0.50 
0.87 


0.32 
0.61 


0.30 
0.74 


0.33 
0.61 


2.30 
2.57 


2.07 
2.28 


2.10 
2.43 


2.26 
2.47 


1.00 
1.00 


0.99 
1.00 


1.00 
1.00 


0.97 
1.00 
















V 


V 


V 


V 


0.90 


0.60 


0.56 


0.61 


2.47 


2.19 


2.19 


2.37 


1.00 


1.00 


1.00 


0.99 
















V 


V 


V 


V 


0.21 


0.12 


0.10 


0.12 


1.87 


2.00 


2.10 


2.07 


0.94 


0.99 


1.00 


1.00 


Mpp A A 'XT 
In \j \j 440 / 


D.U 


Q S 

y .0 


82 4 




-21 47 


±oy .0 


V 
V 


V 
V 


V 
V 


V 
V 


0.15 
0.44 


0.12 
0.25 


0.09 
0.19 


0.09 
0.24 


1.87 
1.92 


1.95 
2.05 


2.07 
2.17 


2.12 
2.12 


0.94 
0.94 


0.94 
1.00 


1.00 
1.00 


1.00 
1.00 
















V 


V 


V 


V 


0.35 


0.24 


0.18 


0.19 


1.90 


2.00 


2.12 


2.14 


0.94 


0.95 


1.00 


1.00 
















V 


V 


V 


V 


0.68 


0.06 


0.06 


0.19 


2.26 


2.32 


2.43 


2.19 


1.00 


1.00 


1.00 


1.00 


NGC 4565 


3.2 


13.3 


135.1 


1019 


-22.51 


244.9 


V 
V 


V 
V 


V 
V 


V 
V 


0.70 
1.20 


0.06 
0.12 


0.06 
0.10 


0.15 
0.36 


2.14 
2.35 


2.26 
2.37 


2.37 
2.55 


2.19 
2.28 


1.00 
1.00 


1.00 
1.00 


1.00 
1.00 


1.00 
1.00 
















V 


V 


V 


V 


1.24 


0.10 


0.10 


0.30 


2.26 


2.37 


2.47 


2.26 


1.00 


1.00 


1.00 


1.00 
















V 


V 


V 


V 


0.29 


0.33 


0.25 


0.30 


2.45 


2.32 


2.24 


2.59 


1.00 


1.00 


1.00 


1.00 


NGC 5470 


3.1 


14.7 


62.2 


158 


-18.01 


109.4 


V 
V 


V 
V 


V 
V 


V 
V 


0.27 
0.53 


0.32 
0.60 


0.16 
0.64 


0.30 
0.56 


2.35 
2.57 


2.21 
2.45 


2.17 
2.43 


2.49 
2.76 


0.97 
1.00 


1.00 
1.00 


low 
l.oq 


1.00 
1.00 
















V 


V 


V 


V 


0.53 


0.59 


0.42 


0.56 


2.49 


2.32 


2.24 


2.63 


1.00 


1.00 


1.00 


1.00 
















V 


V 


V 


V 


0.12 


0.08 


0.15 


0.53 


2.12 


2.21 


2.05 


2.12 


1.00 


1.00 


0.99 


1.00 


in Kj o oyoi 


A Q 


29 2 


loy.o 


1D0 


on x.a 

-ZU.04 


251 1 


V 
V 


V 
V 


V 
V 


V 
V 


0.10 
0.27 


0.08 
0.21 


0.15 
0.27 


0.54 
0.96 


2.10 
2.14 


2.14 
2.14 


2.00 
2.12 


2.05 
2.26 


0.99 
1.00 


1.00 
0.99 


0.96 
1.00 


1.00 
1.00 
















V 


V 


V 


V 


0.23 


0.13 


0.27 


0.98 


2.12 


2.24 


2.07 


2.17 


1.00 


0.99 


0.96 


1.00 
















V 


_L 


V 


V 






0.53 








2.49 








1.00 




PGC 012439 


6.2 


39.4 


0.3 


102 


-19.04 


106.1 


V 
V 
V 


_L 

_L 
± 


V 
V 
V 


V 
V 
V 






0.54 
0.93 
0.96 








2.39 
2.61 
2.53 








1.00 
1.00 
1.00 


















V 


V 


V 


V 




0.61 


0.50 






2.30 


2.32 






1.00 


1.00 




PGC 012798 


7.4 


22.6 


9.7 


161 


-18.97 


101.1 


V 
V 
V 


V 
V 
V 


V 
V 
V 


V 
V 
V 




0.63 
1.10 
1.12 


0.51 
0.90 
0.91 






2.21 
2.39 
2.32 


2.24 
2.45 
2.37 






1.00 
1.00 
1.00 


1.00 
1.00 
1.00 


















V 


V 


V 


V 


0.83 


0.13 


0.13 


0.33 


1.92 


2.07 


2.05 


1.92 


1.00 


1.00 


1.00 


1.00 


PGC 013646 


5.0 


32.6 


34.4 


203 


-19.97 


168.1 


V 
V 


V 
V 


V 
V 


V 
V 


0.84 
1.47 


0.12 
0.24 


0.10 
0.24 


0.30 
0.63 


1.84 
2.00 


2.03 
2.14 


2.05 
2.14 


1.84 
2.00 


1.00 
1.00 


1.00 
1.00 


1.00 
1.00 


1.00 
1.00 
















V 


V 


V 


V 


1.50 


0.25 


0.24 


0.57 


1.92 


2.07 


2.07 


1.92 


1.00 


1.00 


1.00 


1.00 



Table 2 — Continued 



ID 


T 


D 


PA 


-D 25 


B 






Status 




PTo/ 


Ct/ <7% 


fo 








(Mpc) 


(°) 


(") 




(kins" 1 ) 










Pta 












(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) (13) 


(14) (15) 


(16) (17) 


(18) (19) 


(20) (21) 


(22) (23) 
















V 


V 


V 


V 


0.42 


0.68 


1.92 


2.05 


0.99 


1.00 


UGC 09977 


5.3 


31.0 


77.6 


233 


-19.79 


112.7 


V 
V 


V 
V 


V 
V 


V 
V 


0.42 
0.78 


0.69 
1.22 


1.87 
2.00 


1.97 
2.14 


0.96 
1.00 


1.00 
1.00 
















V 


V 


V 


V 


0.78 


1.24 


1.95 


2.07 


0.97 


1.00 
















V 


V 


V 


V 


0.25 


0.18 


2.28 


2.39 


1.00 


1.00 


UGC 10288 


5.3 


32.4 


90.4 


287 


-20.32 


166.6 


V 
V 


V 
V 


V 

± 


V 
V 


0.25 
0.46 


0.18 


2.21 
2.39 


2.32 


1.00 
1.00 


1.00 
















V 


V 


V 


V 


0.45 


0.33 


2.30 


2.43 


1.00 


1.00 
















V 


V 


V 


V 




0.59 




1.73 




1.00 


UGC 10297 


5.1 


39.2 


2.9 


128 


-19.32 


102.8 


< < < 


V 
V 
V 


V 
V 
V 


< < < 




0.46 
1.07 
0.89 




1.64 
1.79 
1.70 




1.00 
1.00 
1.00 



Note. — Galaxy ID (col. 1), morphological type (HypcrLEDA; col. 2), distance (using NASA/IPAC Extragalactic Database mean value of redshift-independent distances except for IC 1553, 
NGC 5470 and PGC 012349 for which wc estimated the distance using HyperLEDA's heliocentric radial velocity from radio measurement and using a Hubble constant Hq — 70 kms -1 ; 
col. 3), PA (col. 4), -D25 (HyperLEDA; col. 5), absolute blue magnitude (HypcrLEDA; col. 6), apparent maximum rotation velocity of gas (HypcrLEDA; col. 7), luminosity profile fij^status 
for projected galactoccntric distances — 0.8r25 < R < — 0.5r25, — 0.5r25 < R < — 0.2r25, 0.2r25 < R < 0.5r25 and 0.5r25 < R < 0.8r25 (V meaning that the fit has been successful, * meaning 
that a bright star made the fit impossible, V meaning that the mid-plane surface brightness was fi > 21.5 mag arcscc -2 thus preventing Am > 4.5, _l_ meaning that although the mid-plane 
brightness was fi < 21.5 mag arcscc -2 , Am < 4.5 and fa meaning that dust lanes not located in the mid-plane arc affecting the fit; cols. 8 to 11), fitted thick to thin disk mid-plane stellar 
density (cols. 12 to 15), fitted thick to thin disk velocity dispersion ratio in the direction of the z axis (cols. 16 to 19) and fitted fraction of the mid-plane light not absorbed by dust (cols. 20 
to 23). For each galaxy, the first row corresponds to a fit using T-r/Yt — 1.2 and without-gas, the second one to Y<r/T t — 1.2 and with-gas, the third one to Tx/Tt — 2.4 and without-gas 
and the fourth one to Y<r/Y t — 2.4 and with-gas. Data in italics indicate fitted values which arc compatible with those of a single disk structure. 



Table 3. Galaxy and fit data (part two) 



ID 




z f 








z T 








St 


/St 






Hi 




A 


m 








(pc 








(pc 


) 












(mag arcscc ) 




(mag arcscc ) 




(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 


(13) 


(14) (15) 


(16) 


(17) 


(18) (19) 


(20) 


(21) 






200 


180 






850 


860 






0.95 


0.99 




24.27 


23.19 




5.5 


4.5 




ESO 079-003 




210 


200 






850 


910 






0.95 


0.92 




24.27 


23.19 




5.5 


4.5 






180 


170 






820 


840 






1.83 


1.87 




24.27 


23.19 




5.5 


4.5 








200 


180 






830 


870 






1.83 


1.80 




24.27 


23.19 




5.5 


4.5 










180 


220 






870 


920 






1.63 


1.69 




25.21 


25.28 






4.5 








190 
160 


230 
220 






880 
820 


930 
900 






1.67 
3.26 


1.70 
3.11 




25.21 
24.71 


25.28 
25.28 




5.5 
5.0 


4.5 
4.5 








180 


220 






860 


900 






3.02 


3.17 




25.21 


25.28 




5.5 


4.5 






190 








850 








0.47 






25.07 






5.0 






F,SD ^44 (197 




210 
190 
210 








900 
840 
890 








0.44 
0.87 
0.81 






25.07 
25.07 
25.07 






5.0 
5.0 
5.0 












130 








340 








2.50 






25.67 






4.5 










140 
120 
130 








340 
340 
340 








2.40 
4.39 
4.29 






25.61 
25.61 
25.61 






A 5 
I 5 
I 5 








180 


160 






690 


590 






1.11 


0.98 




25.65 


24.91 




5.0 


5.0 




IC 0217 




200 


210 






700 


720 






1.09 


0.54 




25.65 


25.41 




5.0 


5.5 






170 


140 






670 


580 






2.12 


1.86 




25.65 


24.91 




5.0 


5.0 








180 


200 






680 


720 






2.09 


1.02 




25.65 


25.41 




5.0 


5.5 








170 


120 






680 


610 






1.56 


2.22 




25.45 


25.42 




4.5 


4.5 




IC 1197 




180 


130 






690 


610 






1.54 


2.34 




25.45 


25.42 




4.5 


4.5 






160 


110 






670 


600 






2.89 


3.98 




25.45 


25.42 




4.5 


4.5 








170 


120 






670 


600 






2.90 


4.12 




25.45 


25.42 




4.5 


4.5 








200 


220 


240 




680 


910 


1280 




1.03 


0.33 


0.36 


24.23 


24.49 


25.59 


4.5 


5.5 


5.5 


IC 1553 




310 


240 


250 




1320 


1040 


1280 




0.27 


0.27 


0.37 


25.23 


25.49 


25.59 


5.5 


6.5 


5.5 




190 


210 


230 




660 


870 


1250 




1.92 


0.68 


0.69 


24.23 


24.49 


25.59 


4.5 


5.5 


5.5 






310 


230 


250 




1400 


1020 


1270 




0.48 


0.55 


0.70 


24.73 


25.49 


25.59 


5.0 


6.5 


5.5 






190 


180 






990 


1030 






0.50 


0.43 




24.49 


25.00 




5.0 


5.5 




IC 1970 




220 


200 






1070 


1050 






0.47 


0.44 




24.49 


25.00 




5.0 


5.5 






190 


180 






990 


1040 






0.93 


0.78 




24.49 


25.00 




5.0 


5.5 








200 


200 






1030 


1050 






0.91 


0.79 




24.49 


25.00 




5.0 


5.5 








160 


120 






490 


490 






1.39 


1.87 




25.59 


25.73 




4.5 


4.5 




IC 2058 




170 


120 






500 


490 






1.36 


1.97 




25.59 


25.73 




4.5 


4.5 






150 


110 






490 


480 






2.67 


3.43 




25.59 


25.73 




4.5 


4.5 








160 


110 






490 


480 






2.66 


3.56 




25.59 


25.73 




4.5 


4.5 





Table 3 — Continued 



ID 






Zt 






z T 










/St 






W 








A 


m 








(PC) 






(PC) 














(mag arcscc 2 ) 






(mag arcscc 2 ) 




(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 


(13) 


(14) 


(15) 


(16) 


(17) 


(18) 


(19) 


(20) 


(21) 






230 


220 


370 




930 


840 


1700 




1.42 


1.35 


0.60 




24.65 


24.86 


25.92 




4.5 


5.0 


5.0 


IC 2135 


270 


250 


230 


390 


990 


940 


850 


1720 


1.56 


1.44 


1.37 


0.62 


25.35 


24.65 


24.86 


25.92 


4.5 


4.5 


5.0 


5.0 




220 


210 


340 




910 


830 


1640 




2.67 


2.53 


1.20 




24.65 


24.86 


25.92 




4.5 


5.0 


5.0 






230 


230 


370 




910 


840 


1680 




2.68 


2.54 


1.16 




24.65 


24.86 


25.92 




4.5 


5.0 


5.0 




170 


140 






530 


470 






2.10 


2.47 






26.00 


25.53 






4.5 


4.5 






IC 5052 


190 


140 






530 


470 






2.05 


2.57 






26.00 


25.53 






4.5 


4.5 






160 
1 70 


130 
130 






530 
530 


470 
480 






3.75 
3.86 


4.42 
4-59 






26.00 
26.00 


25.53 
25.53 






4.5 
4.5 


4.5 
4.5 








230 


190 


160 


220 


930 


810 


730 


850 


0.78 


0.82 


1.00 


0.85 


25.08 


24.75 


24.13 


25.51 


5.0 


5.5 


5.0 


5.5 


NGC 0522 


260 


210 


190 


240 


980 


830 


760 


880 


0.69 


0.77 


0.92 


0.82 


25.08 


24.75 


24.13 


25.51 


5.0 


5.5 


5.0 


5.5 


190 


160 


150 


210 


800 


740 


710 


840 


1.98 


1.83 


1.85 


1.65 


24.58 


24.25 


24.13 


25.51 


4.5 


5.0 


5.0 


5.5 




240 


200 


170 


220 


950 


830 


740 


850 


1.37 


1.39 


1.77 


1.59 


25.08 


24.75 


24.13 


25.51 


5.0 


5.5 


5.0 


5.5 






450 








1760 








0.69 








25.18 








4.5 






NGC 0678 




580 
420 
540 








2890 
1710 
2640 








0.37 
1.32 
0.77 








25.68 
25.18 
25.68 








5.0 
4.5 
5.0 










120 


110 






580 


530 






1.49 


1.55 






24.90 


25.48 






4.5 


5.0 




NGC 1351A 




130 


120 






590 


530 






1.47 


1.58 






24.90 


25.48 






4.5 


5.0 






120 
130 


110 
110 






570 
570 


510 
520 






2.70 
2.79 


2.83 
2.91 






24.90 
24.90 


25.48 
25.49 






4.5 
4.5 


5.0 
5.0 








100 


120 






370 


350 






1.46 


1.71 






24.61 


24.80 






5.0 


5.0 




NGC 1495 




100 


130 






380 


350 






1.47 


1.70 






24.61 


24.80 






5.0 


5.0 






90 
90 


110 
120 






370 
370 


340 
340 






2.72 
2.72 


3.20 
3.15 






24.61 
24.61 


24.80 
24.80 






5.0 
5.0 


5.0 
5.0 








180 


90 






600 


330 






0.38 


2.95 






25.75 


25.54 






5.0 


4.5 




NGC 1827 




210 


100 






700 


330 






0.30 


3.04 






25.75 


25.54 






5.0 


4.5 






120 
200 


90 
90 






370 
710 


330 
330 






2.69 
0.55 


5.21 
5.39 






25.25 
25.75 


25.54 
25.54 






4.5 
5.0 


4.5 
4.5 








210 


180 






610 


560 






0.95 


1.15 






25.70 


25.64 






5.5 


5.5 




NGC 3501 




240 


200 






650 


560 






0.81 


1.12 






25.70 


25.64 






5.5 


5.5 






210 
220 


170 
180 






600 
630 


550 
560 






1.84 
1.67 


2.16 
2.09 






25.70 
25.70 


25.64 
25.64 






5.5 
5.5 


5.5 
5.5 








290 


300 






1470 


1520 






1.67 


1.73 






24.41 


25.01 






5.0 


5.5 




NGC 3628 




300 


330 






1510 


1540 






1.66 


1.82 






24.41 


25.01 






5.0 


5.5 






270 
280 


290 
310 






1460 
1470 


1510 
1510 






3.05 
3.04 


3.12 
3.26 






24.41 
24.41 


25.01 
25.01 






5.0 
5.0 


5.5 
5.5 





Table 3 — Continued 



ID 




Zt 








z T 








St 


/St 






W 








A 


m 








(PC) 






(PC) 














(mag arcscc 2 ) 






(mag arcscc 2 ) 




(1) 


(2) 


(31 


(41 


(5) 


(6) 


(71 


(81 


(9) 


(10) 


(ii) 


(12) 


(13) 


(14) 


(1 51 


(161 


(17) 


(18) 


(191 


(201 


(21) 






120 


110 


180 




720 


680 


1190 




1.12 


1.20 


1.09 




22.99 


22.93 


24.67 




4.5 


4.5 


4.5 


NGC 4013 




130 


120 


190 




730 


700 


1210 




1.11 


1.20 


1.08 




22.99 


22.93 


24.67 




4.5 


4.5 


4.5 






120 


110 


180 




700 


680 


1140 




2.19 


2.26 


2.22 




22.99 


22.93 


24.67 




4.5 


4.5 


4.5 




150 


190 


170 


190 


700 


730 


710 


820 


1.47 


0.84 


0.87 


0.98 


25.54 


25.65 


25.44 


25.93 


5.0 


5.5 


5.5 


5.5 


NGC 4330 


170 


200 


190 


200 


710 


760 


740 


830 


1.48 


0.84 


0.82 


0.99 


25.54 


25.65 


25.44 


25.93 


5.0 


5.5 


5.5 


5.5 


150 


180 


140 


170 


680 


710 


620 


800 


2.74 


1.71 


2.19 


1.89 


25.54 


25.65 


24.94 


25.93 


5.0 


5.5 


5.0 


5.5 




150 


190 


180 


190 


690 


730 


710 


820 


2.81 


1.66 


1.55 


1.88 


25.54 


25.65 


25.44 


25.93 


5.0 


5.5 


5.5 


5.5 




330 


300 


300 


350 


1040 


1110 


1200 


1400 


0.47 


0.30 


0.28 


0.32 


25.98 


25.97 


25.98 


25.87 


5.0 


6.0 


6.0 


5.0 


NGC 4437 


380 


320 


330 


390 


1200 


1120 


1300 


1630 


0.36 


0.31 


0.25 


0.26 


25.98 


25.97 


25.98 


25.87 


5.0 


6.0 


6.0 


5.0 


310 


280 


290 


350 


990 


1050 


1190 


1330 


0.99 


0.65 


0.54 


0.64 


25.98 


25.97 


25.98 


25.87 


5.0 


6.0 


6.0 


5.0 




340 


300 


320 


370 


1100 


1080 


1250 


1500 


0.81 


0.62 


0.51 


0.56 


25.98 


25.97 


25.98 


25.87 


5.0 


6.0 


6.0 


5.0 




240 


420 


450 


380 


930 


2080 


2290 


1600 


1.85 


0.19 


0.21 


0.55 


25.94 


25.70 


25.74 


25.84 


4.5 


6.0 


6.0 


4.5 


NGC 4565 


250 


460 


460 


430 


950 


2140 


2350 


1850 


1 88 


20 


21 


45 


25 94 


25.70 


25.74 


25 84 


4 5 


6.0 


6.0 


4.5 


220 


410 


430 


360 


930 


2010 


2360 


1560 


3.33 


0.39 


0.38 


1.04 


25.94 


25.70 


25.74 


25.84 


4.5 


6.0 


6.0 


4.5 




240 


440 


450 


410 


930 


2200 


2400 


1740 


3.40 


0.36 


0.39 


0.91 


25.94 


25.70 


25.74 


25.84 


4.5 


6.0 


6.0 


4.5 




90 


70 


90 


90 


430 


330 


380 


460 


0.92 


0.98 


0.73 


1.03 


25.58 


24.64 


24.70 


25.99 


5.0 


5.0 


5.0 


5.0 


NGC 5470 


100 


90 


110 


100 


440 


350 


460 


470 


0.87 


0.92 


0.49 


1.04 


25.58 


24.64 


25.70 


25.99 


5.0 


5.0 


6.0 


5.0 


80 


70 


70 


80 


420 


330 


310 


440 


1.72 


1.83 


1.94 


1.97 


25.58 


24.64 


24.20 


25.99 


5.0 


5.0 


4.5 


5.0 




90 


80 


90 


90 


420 


330 


390 


450 


1.73 


1.75 


1.23 


1.95 


25.58 


24.64 


24.70 


25.99 


5.0 


5.0 


5.0 


5.0 




310 


270 


240 


190 


1230 


1190 


880 


700 


0.33 


0.22 


0.39 


1.35 


25.80 


25.87 


25.93 


25.70 


5.0 


6.5 


6.5 


5.0 


NGC 5981 


330 


280 


250 


200 


1310 


1200 


890 


710 


0.30 


0.23 


0.40 


1.38 


25.80 


25.87 


25.93 


25.70 


5.0 


6.5 


6.5 


5.0 


290 


250 


220 


170 


1110 


970 


880 


690 


0.73 


0.57 


0.72 


2.58 


25.80 


25.37 


25.93 


25.70 


5.0 


6.0 


6.5 


5.0 




310 


280 


240 


180 


1220 


1210 


890 


700 


0.63 


0.42 


0.73 


2.58 


25.80 


25.87 


25.93 


25.70 


5.0 


6.5 


6.5 


5.0 








130 








620 








1.65 








25.41 








4.5 




PGC 012439 






110 








640 








1.68 








25.41 








4.5 








































120 








620 








2.97 








25.41 








4.5 










130 








620 








3.07 








25.41 








4.5 








130 


130 






540 


560 






1.73 


1.43 






25.70 


25.66 






4.5 


4.5 




PGC 012798 




140 


140 






540 


570 






1.76 


1.47 






25.70 


25.66 






4.5 


4.5 






120 


120 






530 


550 






3.12 


2.67 






25.70 


25.66 






4.5 


4.5 








130 


130 






540 


550 






3.20 


2.71 






25.70 


25.66 






4.5 


4.5 






180 


250 


240 


230 


570 


960 


950 


750 


1.84 


0.36 


0.35 


0.76 


25.40 


24.93 


25.44 


26.00 


4.5 


5.5 


6.0 


5.5 


PGC 013646 


200 


270 


280 


260 


580 


1010 


1070 


790 


1.83 


0.32 


0.29 


0.69 


25.40 


24.93 


25.94 


26.00 


4.5 


5.5 


6.5 


5.5 


180 


250 


250 


220 


570 


970 


970 


740 


3.36 


0.65 


0.65 


1.49 


25.40 


24.93 


25.44 


26.00 


4.5 


5.5 


6.0 


5.5 




190 


250 


260 


240 


580 


950 


980 


770 


3.35 


0.69 


0.65 


1.34 


25.40 


24.93 


25.44 


26.00 


4.5 


5.5 


6.0 


5.5 



Table 3 — Continued 



ID 






Zt 






z T 








/St 




W 




Am 






(PC) 






(PC) 








(mag arcscc 2 ) 


^mag 


arcscc ) 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) (11) 


(12) (13) 


(14) (15) 


(16) (17) 


(18) (19) 


(20) (21) 






280 


230 






870 


800 




0.96 


1.64 


25.79 


25.93 


4.5 


4.5 


UGC 09977 




280 


240 






880 


810 




0.96 


1.66 


25.79 


25.93 


4.5 


4.5 






250 


210 






860 


770 




1.83 


3.03 


25.79 


25.93 


4.5 


4.5 






280 


230 






860 


790 




1.83 


3.07 


25.79 


25.93 


4.5 


4.5 






230 


270 






1030 


1310 




0.75 


0.57 


25.58 


25.71 


5.5 


5.5 


UGC 10288 




260 


300 






1050 


1340 




0.76 


0.59 


25.58 


25.71 


5.5 


5.5 






230 








1010 






1.40 




25.58 




5.5 








240 


280 






1040 


1300 




1.36 


1.10 


25.58 


25.71 


5.5 


5.5 ■ • • 








250 








680 






1.16 




25.89 




4.5 


UGC 10297 






290 








720 






0.89 




25.89 




4.5 


























240 








670 






2.16 




25.89 




4.5 








280 








720 






1.74 




25.89 




4.5 



Note. — Galaxy ID (col. 1), fitted thin disk scalchcight for projected galactocentric distances — 0.8r25 < R < — 0.5r25, — 0.5r25 < R < — 0.2r25, 0.2r25 < R < 0.5r25 and 
0.5r25 < R < 0.8r25 (cols. 2 to 5), fitted thick disk scalchcight (cols. 6 to 9), fitted thick to thin stellar column mass density (cols. 10 to 13), limiting magnitude of the fit 
(cols. 14 to 17) and dynamical range of the fit (cols. 18 to 21). For each galaxy, the first row corresponds to a fit using Tx/Tt — 1.2 and without-gas, the second one to 
Tt/Ti — 1.2 and with-gas, the third one to T-r/Yt — 2.4 and without-gas and the fourth one to Tx/T t — 2.4 and with-gas. Data in italics indicate fitted values which are 
compatible with those of a single disk structure. 



